ESTI  FILE  COPY 


Spin-Wave  Effects 
in  the  Magnetization  Reversal 
of  a  Thin  Ferromagnetic  Film 


27  August  1964 


vndar  Electronic  Syeteme  Division  Contrsct  AF  19  (628)-500  by 


Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lexington,  Mnnnnchunettn 


>-'v 

■,vv%^.v  -> 

%,'.lf-.*r‘  *  •■'  «  i 


■i 


The  work  reported  in  this  document  was  performed  at  Lincoln 
Laboratory,  a  center  for  research  operated  by  Massachusetts 
Institute  of  Technology,  with  the  support  of  the  U.S.  Air  Force 
under  Contract  AF  19  (628)- 500. 


Non-Lincoln  Recipients 

PLEASE  DO  NOT  RETURN 


Permission  is  given  to  destroy  this  document 
when  it  is  no  longer  needed. 


tA'.-  v- 

.  •  u-f'i.-.  vlr  .-  -tv* 


■-■V' 


90 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


SPIN-WAVE  EFFECTS  IN  THE  MAGNETIZATION  REVERSAL 
OF  A  THIN  FERROMAGNETIC  FILM 


K.  J.  HARTE 
Group  24 


TECHNICAL  REPORT  364 

27  AUGUST  1964 


LEXINGTON 


MASSACHUSETTS 


SPIN-WAVE  EFFECTS  IN  THE  MAGNETIZATION  REVERSAL 
OF  A  THIN  FERROMAGNETIC  FILM* 


ABSTRACT 


The  influence  of  spin  woves  on  rapid  ratotionol  mognetizotion  reversol  (switching)  in  o  thin  ferromognetic 

film  is  investigated  by  meons  of  o  semiclossicol ,  continuum  theory  which  includes  externol,  onisotropy, 

exchange,  ond  mognetostotic  (dipolor)  fields.  To  simplify  the  mognetostotic  field,  o  "thin-film  opproxi- 

motion"  is  introduced,  in  which  the  mognetizotion  is  reploced  by  its  overoge  over  the  film  thickness. 

From  o  stochostic  model  for  the  microstructure  of  o  polycrystoll ine  film,  the  equilibrium  mognetizotion 

configurotion  M(T)is  derived.  Plonor  fluctuotions  of  M  from  its  meon  direction  ore  found  which 

o 

hove  the  c  ha  roc  ter  i  sties  of  “longitudinol  ripple,"  nomely,  wove  vectors  k  porollel  to  m  ond  wovelengths 

-4  .  .  .  °  . 

greoter  thon  on  exchange  cutoff  2itX  ^  10  cm.  Campanents  with  wove  vectors  in  directions  other  thon 
±m^  ore  ottenuoted  by  mognetostotic  forces,  while  exchange  forces  ottenuote  components  with  wove¬ 
lengths  less  than  The  magnetizotion  dispersion  6  [  rms  ongulor  deviotion  of  M(T)  from  m^]  is  olso 

colculoted.  A  brief  discussion  is  given  of  the  uniform  rototionol  switching  mode  (without  spin  woves), 
with  particulor  ottention  to  undamped  ond  overdamped  cases.  From  the  spin-wove  equotions  of  motion, 
the  spectrum  applicable  to  a  parollel  resononce  situotion  (externol  field  in  the  film  plone)  is  first  ob- 
toined,  ond  long-wovelength  mognetostotic  distortion  is  noted.  Then  the  tronsient  spin-wove  response  is 
colculoted  for  o  switching  field  H  suddenly  turned  on  ot  t  =  0.  It  is  found  that  if  m  (t)  rototes  foster  than 

1  ^  P  ^  o 

longitudinal  spin  woves  [k  ||  m  (0)]  con  relox,  the  mognetizotion  goes  through  a  transient  stote  of  high 

.  .  2  . 

mognetostotic  energy,  ond  o  spin-wove  reaction  torque  (proport iorra I  to  5  )  is  exerted  on  the  uniform 

mode.  If  H^  is  less  rhan  o  critical  field  reaction  torque  at  some  point  in  the  switching  process 

is  greoter  than  the  reversing  torque  and  the  uniform  mode  becomes  locked;  rototionol  reversal  connot 

proceed  until  initiolly  longitudinol  spin  woves  hove  relaxed  into  components  propogoting  in  the  inston- 

toneous  direction  of  ni^(0*  Such  o  highly  damped  process  is  suggestive  of  the  noncoherent  reversol  mode 

observed  in  thin  films.  For  H  >  H  ,  reversal  takes  place  by  o  modified  uniform  rototion;  H  may 

P  pc  ^  .  pc 

therefore  be  identified  as  the  threshold  field  for  coherent  rototion.  The  colculoted  dependence  of  H  on 

pc 

o  bios  field  comperes  fovorobly  with  experiment.  The  S-dependence,  if  8  con  be  meosured  independ¬ 
ently,  should  provide  o  cruciol  test  of  the  theory. 
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SPIN-WAVE  EFFECTS  IN  THE  MAGNETIZATION  REVERSAL 
OF  A  THIN  FERROMAGNETIC  FILM 


I.  INTRODUCTION 

A.  MAGNETIZATION  REVERSAL  AND  THIN  FILMS: 

EXPERIMENTAL  BACKGROUND 

In  the  past  decade  considerable  interest  has  arisen  in  the  phenomenon  of  magnetization 
reversal  in  ferromagnetic  media,  i.e.,  an  irreversible  change  in  the  direction  of  magnetization 
following  the  application  of  a  pulse  magnetic  field.  Although  this  interest  received  its  initial 
impetus  from  the  purely  practical  possibility  of  ferromagnetic  storage  elements  for  high-speed 
computers,  it  soon  became  apparent  that  the  physical  processes  involved  were  either  imperfectly 
understood  or  else  completely  unknown.  The  basic  distinction  between  magnetization  reversal, 
or  switching,  on  the  one  hand,  and  various  magnetic  resonance  phenomena,  which  are  generally 
much  better  understood,  on  the  other,  is  that  the  former  involves  large  displacements  of  the  mag¬ 
netization  M  from  its  equilibrium  direction.  How  M  returns  to  equilibrium  is  the  problem  to 
which  we  shall  devote  ourselves  in  this  work. 

In  particular,  we  shall  be  concerned  with  magnetization  reversal  in  thin  films,  for  the  re¬ 
sultant  theoretical  simplicity  seems  well  worth  the  loss  in  generality.  Furthermore,  experimental 
results  with  thin  films  tend  to  be  consistent  and  unambiguous,  in  contrast  to  the  somewhat  confused 
situation  with  respect  to  switching  in  bulk  ferromagnets .  Among  the  significant  properties  of 
thin  films  are:  (a)  the  shape  demagnetizing  field  confines  M  to  the  film  plane  or  nearly  so  (in 
the  absence  of  a  strong  field  or  easy  axis  normal  to  the  plane);  (b)  eddy-current  effects  in  ferro¬ 
magnetic  metals  may  be  eliminated  by  choosing  a  film  thickness  much  less  than  the  skin  depth; 

(c)  the  equilibrium  magnetic  configuration  may  be  directly  observed  since  it  is  invariably  the 

o 

same  as  the  surface  configuration;  (d)  the  typically  small  crystallite  size  (~300  A)  of  polycrys- 
tallinc  films,  in  conjunction  with  (a),  results  in  a  remanent  state  which  is  very  nearly  single¬ 
domain;  (e)  single-crystal  films  of  various  ferromagnetic  metals  and  nonmetals  may  be  studied 
(these  are  usually  epitaxially  grown);  (f)  low  anisotropy  (a  few  oersteds)  is  readily  obtained, 
thus  allowing  magnetization  reversal  to  take  place  in  low,  fast-rise  fields. 

Since  the  early  experiments  of  Blois,  many  workers  have  studied  pulse  switching  in  thin 
2-10 

films,  usually  of  Permalloy  (Ni-Fe  alloy  near  the  zero  magnetostriction  0.83-0.17  com¬ 

position),  and  possessing  uniaxial  anisotropy  with  easy  axis  in  the  film  plane.  In  very  idealized 
form,  these  experiments  have  been  performed  in  the  following  manner  (see  the  above  references 
for  details,  especially  Refs.  5  and  9).  A  fast  rise-time  (^10  sec)  pulse  field  of  a  few  oer¬ 
steds  magnitude  is  applied  in  the  film  plane  at  an  obtuse  angle  to  the  mean  magnetization  m  of 

t  ^ 

an  initially  single-domain  film.  The  switching  time  t^  and  other  information  on  the  reversal 

fMany  of  the  early  experiments  were  performed  with  Hp  antiparallel  ta  mQ,  resulting  often  in  bi-directianal  re¬ 
versal,  with  its  awn  peculiar  characteristics.  We  shall  exclude  this  limiting  case  from  aur  treatment. 
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PULSE  SWITCHING  FIELD  (Hp) 


Fig.  1.  Schemotic  illustrotion  of  o  thin-film 
switching  curve  showing  (1 )  low-speed  region, 
(2)  intermediote-speed  region,  ond  (3)  high¬ 
speed  region.  Also  shown  ore  threshold  fields 
for  eoch  region. 


observed  by  Olson  and  Pohm  ).  The  reversal  process  in  region  1  (low  speed;  ~  10  ^  sec) 


process  is  obtained  from  the  voltage  induced  in  a  pickup  loop  with  axis  parallel  to  11^.  This 
voltage  is  termed  the  longitudinal  switching  signal;  further  information  may  be  obtained  from 
the  transverse  switching  signal,  induced  in  an  orthogonal  pickup  loop  with  axis  also  in  the  film 
plane 

The  main  results  of  these  experiments  are  conveniently  displayed  in  the  form  of  switching 

curves,  in  which  inverse  switching  time  is  plotted  as  a  function  of  the  pulse  field.  Figure  1 

shows  a  typical  thin-film  switching  curve;  we  note  the  presence  of  three  distinct  regions  (first 

5,  - -  - -  ^  >  .„-6 

s 

has  been  identified  as  the  nucleation  of  reverse  domains  and  subsequent  domain-wall  motion. 

The  threshold  field  for  this  process,  H  .,  is  the  wall  coercive  force;  a  fairly  successful  model 

11 

for  domain-wall  switching  has  been  proposed  by  Conger  and  Essig.  We  shall  not  discuss 

region  1  further  since  the  physical  processes  involved  are  well  understood. 

-8 

We  proceed  now  to  region  3  (high  speed:  t^  <  10  sec),  where  the  reversal  process  is 
clearly  a  coherent,  or  nearly  coherent,  rotation  of  M.  A  phenomenological  theory  of  such  a 

process,  based  on  a  damped  gyromagnetic  equation  such  as  that  proposed  by  Landau  and 

12  4  5 

Lifshitz,  has  been  found  by  Smith,  and  by  Olson  and  Pohm,  to  describe  switching  in  region  3 

fairly  well.  However,  there  are  two  difficulties.  First,  and  most  disturbing,  the  threshold 

field  is  always  found  to  be  greater  (by  a  factor  ~1.5)  than  the  theoretical  Stoner-Wohlfarth 

threshold  H  for  a  uniformly  magnetized  uniaxial  film.  Second,  the  damping  of  high-speed 
pi 

reversal  is  found  to  be  3  or  4  times  the  intrinsic  damping,  as  obtained  from  low-power  ferro- 

10  14  15  10  14 

magnetic  resonance  linewidth  "  "  or  decay  of  free  oscillations  of  the  magnetization. 

1 5 

Smith  and  Harte  have  suggested  that  the  discrepancy  in  the  threshold  field  may  be  caused 
by  dispersion  in  the  direction  and  amplitude  of  the  uniaxial  anisotropy . ^ ^  By  assuming  that 
those  regions  of  a  film  with  the  highest  threshold  control  coherent  rotation,  they  have  been  able 
to  account  for  observed  thresholds  with  fairly  plausible  values  of  anisotropy  dispersion.  However, 
their  assumption  is  difficult  to  justify,  since  the  effective  anisotropy  in  coherent  rotation  of  M 
should  be  the  mean,  or  close  to  it.  Furthermore,  interactions  between  regions  were  neglected. 
The  discrepancy  in  damping,  for  which  no  explanation  exists,  implies  that  there  is  a  loss 
mechanism  operative  in  (large-angle)  switching  that  is  not  present  in  (small-angle)  resonance. 


13 


tThe  exoct  meoning  of  "switching  time"  vories  somewhot  from  worker  to  worker  ond  moy  depend  upon  the  noture 
of  the  reversal  process.  For  present  purposes,  t^  need  not  be  precisely  defined  and  moy  just  be  considered  os 
some  meosure  of  the  durotion  of  the  longitudinol  switching  signol . 


2 


It  should  now  be  clear  that  even  such  a  conceptually  simple  process  as  coherent  domain  rotation 
is  only  imperfectly  understood. 

Region  2  (intermediate  speed:  10  t^  £  10  )  has  baffled  workers  since  its  discovery  * 

and  the  physical  processes  involved  are  not  known.  The  threshold  field  H  however,  has  been 

4,5 


found  to  coincide  with  the  theoretical  rotational  threshold  H 


Pt’ 


This,  in  addition  to  the  obser¬ 


vation  of  a  sizable  transverse  switching  signal,  implies  that  reversal  takes  place  at  least  partially 
by  a  rotational  mechanism.  However,  switching  in  region  2  is  about  twenty  times  slower  than 
that  predicted  by  the  above-mentioned  theory  of  coherent  reversal.  Humphrey^  has  suggested 
that  reversal  starts  by  simple  coherent  rotation  but  quickly  breaks  up  into  what  he  has  termed 
noncoherent  rotation.  Although  this  picture  is  consistent  with  all  experimental  evidence,  in¬ 
cluding  the  shape  of  the  longitudinal  switching  signal  (an  initial  spike  followed  by  a  long  tail^), 
a  physical  model  was  lacking,  and  it  was  not  even  clear  what  precipitated  the  breakup. 

g 

Humphrey  and  Gyorgy  have  obtained  a  phenomenological  description  of  noncoherent  roti^cion 
1 8 

based  on  the  Gilbert  modification  of  the  Landau-Lifshitz  equation  with  damping  about  100  times 

1 9 

the  intrinsic  damping,  but  no  explanation  of  this  extraordinarily  large  loss  was  given.  Harte  has 
proposed  a  model  of  noncoherent  rotation  based  on  angular  anisotropy  dispersion,  but  magneto¬ 
static  interactions  were  neglected  and  the  remaining  effects  are  far  too  small  to  account  for  the 
large  loss.  Smith  and  Harte^^  have  suggested  that  intermediate-speed  switching  may  take  place 
by  a  sequential  process,  observed  quasistatically  by  Smith,^^  which  they  termed  labyrinth  propa¬ 
gation.  (An  incipient  labyrinth  may  be  seen  in  the  photograph  in  P"ig.  2  as  a  reverse  domain  which 

has  propagated  from  the  lower  right  film  edge.)  Another  possibility  is  a  partial  rotation  process, 

21 

observed  quasistatically  by  Methfessel,  et  al.  However,  no  attempt  was  made  in  either  case 
to  calculate  the  dynamic  characteristics  of  such  processes,  and  it  is  not  known  how  valid  these 
models  are. 


LONGITUDINAL  RIPPLE 
(V-M  «  0) 


250-A  PERMALLOY  FILM 


^ 

XXXX 

y^  y^  y^  y^ 

TRANSVERSE  RIPPLE 
(7-M  ^  0) 


RANDOM  ANISOTROPY  MODEL 


Fig.  2.  Magnetization  ripple  in  a  thin  film.  Tap  center:  defacused  electran  micragraph 
af  a  250-A  Permollay  film  (courtesy  of  M.S.  Cahen)  showing  longitudinal  ripple  (fine 
structure)  and  large-angle  domain  walls  (heavy  black  ar  white  lines).  Tap  left  and  right; 
schematic  illustration  of  longitudinal  and  transverse  ripple,  respectively.  Bottom:  anisot¬ 
ropy  model  showing  randomly  varying  and  uniform  uniaxial  anisotropies. 
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In  summary,  the  unsolved  problems  in  thin-film  switching  are:  (a)  accounting  for  the  ele¬ 
vated  coherent  rotational  threshold;  (b)  discovering  the  reason  for  breakup  of  coherent  rotation  in 
region  2;  (c)  finding  the  origin  of  the  anomalous  damping  in  coherent  rotation;  and  (d)  uncovering 
the  physical  mechanism  of  noncoherent  rotation.  The  theory  we  develop  provides  answers  to  (a) 
and  (b),  provides  a  qualitative  answer  to  (c)  —  further  calculation  should  result  in  a  quantitative 
comparison  of  theory  and  experiment  on  this  point,  and  leads  to  a  suggestion  for  (d)  which,  although 
plausible,  cannot  be  fairly  judged  without  additional  development  of  the  model. 

B.  THEORETICAL  FOUNDATION 

Our  analysis  of  magnetization  reversal  is  based  on  a  semiclassical,  continuum  theory  in 
which  the  quantum  mechanical  spin-density  operator  of  the  ferromagnetic  electrons  is  replaced 
by  a  constant  times  a  classical  vector  field  M(r,  t)  —  the  magnetic  moment  density  or  magnetiza¬ 
tion.  This  replacement  is  possible  provided  the  wavelengths  of  all  significant  disturbances  of 
the  spin  system  are  much  greater  than  the  lattice  constant  and  provided  wc  are  not  too  near  the 

Curie  temperature.  For  a  thorough  discussion  of  the  validity  of  this  semiclassical  theory,  the 

22 

reader  is  referred  to  a  paper  by  Herring  and  Kittel.  We  assume  that  all  fluctuations  of  M(r,t) 
from  its  spatial  average  rn^(t)  are  small,  i.e.,  that  the  sample  is  nearly  a  single  domain.  This 
condition  we  write  as 

I  6M  I  «  M  (I- la) 

I  I  o 

where 


M(r,  t)  =  ^  6M(r,  t) 


(I-lb) 


and 


M  =  ImI  .  (I-lc) 

O  I  I  '  / 

The  fluctuations  61VI  will  be  assumed  constant  across  the  film  thickness,  i.c.,  with  wave 
vectors  k  in  the  film  plane.  However,  generalization  to  three-dimensional  k-space  is  not 
difficult  so  that  it  may  be  possible  to  generalize  the  theory  to  include  reversal  processes  in 
bulk  ferromagnets. 

1 2 

As  suggested  phenomenologically  by  Landau  and  Lifshitz,  and  later  established  on  a  very 

22 

general  basis  by  Herring  and  Kittel,  the  equation  of  motion  for  M  may  be  written  in  the  form 


=  — yM  X  H^^^  +  damping  term  (1-2) 

where  y  is  the  absolute  value  of  the  gyromagnetic  ratio  and  H^^^(r,  t)  is  an  effective  magnetic 
field  which,  for  our  purposes,  consists  of  an  exchange  field  H^,  a  magnetostatic  field  H^,  an 
effective  anisotropy  field  H^,  and  an  external  field  H^^.  In  addition,  H^^^  includes  the  effects  of 
magnetostriction  and  finite  conductivity.  In  the  theory  we  shall  develop,  damping  processes  are 
of  secondary  importance;  therefore  the  form  of  the  damping  term  in  (1-2)  is  not  crucial  to  our 
work.  In  Sec.  IV  we  shall  briefly  discuss  the  various  damping  terms  which  have  been  proposed, 
but  for  most  of  this  work  damping  will  be  neglected. 

The  constant  y  is  given  by 


y  = 


ge 

Imc 


(1-3) 


4 


where  e  and  m  are  the  electronic  charge  and  mass,  c  is  the  velocity  of  light,  and  the  spectro¬ 
scopic  splitting  factor  g  is  slightly  greater  than  2  due  to  incomplete  quenching  of  orbital  angular 
momentum.^^ 

The  exchange  field  H  ,  which  arises  from  the  short-range  interactions  responsible  for 
^  22 

ferromagnetism,  has  been  shown  by  Herring  and  Kittel  to  be  given,  for  a  cubic  lattice,  by 

H  =  — ^  ^  U  (1-4) 

e  n  /t  ^ 

M 

o 

where  A  is  an  exchange  constant  (in  erg/cm).  It  should  be  noted  that  as  given  by  (1-4),  is 
independent  of  atomic  model  and  follows  from  the  energy  expression  of  lowest  order  in  the  spatial 
derivatives  of  IVl  allowed  by  the  symmetry  of  the  lattice,  and  invariant  with  respect  to  spin  ro¬ 
tations  and  inversions. 

The  magnetostatic  field  11^  has  its  origin  in  dipole-dipole  interactions  and  in  our  continuum 
theory  may  be  found  from  Maxwell's  magnetostatic  equations 

V  .  (H  +  4irM)  =  0  (I-5a) 

m 


V  X  H  =  0 
m 


(I- 5b) 


We  devote  Sec.  II  to  the  solution  of  (I-5a,  b)  with  boundary  conditions  appropriate  to  thin-film 
geometry. 

Anisotropy  forces  in  thin  polycrystalline  films  are  of  two  general  types,  uniform  and  local. 

A  uniform  uniaxial  anisotropy  is  usually  found  to  be  induced  along  the  magnetization  direction 
during  film  deposition.^  Although  the  physical  origins  of  this  anisotropy  are  not  completely  under¬ 
stood  (for  a  discussion  of  these  the  reader  is  referred  to  a  paper  by  Smith  ),  a  satisfactory 
phenomenological  description  may  be  obtained  from  the  energy  expression 

E  =  -K  M  "^(M  •  1  (1-6) 

a  o  o  X 

where  i^  is  a  unit  vector  along  the  easy  axis  and  is  a  positive  constant.  Local  anisotropy 
forces  arise  from  inhomogeneities  and  will  be  discussed  in  the  introduction  to  Sec.  III.  Magneto- 
strictive  effects,  whether  uniform  or  local,  will  be  included  in  the  effective  anisotropy  field. 

External  fields  may  be  steady  or  time-varying  and  entail  no  conceptual  problems.  Rise 
times  of  the  order  of  10  sec  for  a  pulse  field  of  a  few  oersteds  are  now  readily  available, 
and  we  shall  assume  that  pulse  switching  fields  may  be  represented  by  step  functions. 

24 

The  primary  effect  of  finite  conductivity  in  metallic  films  is  eddy-current  damping.  Smith 
has  shown  that  the  mean  eddy-current  contribution  to  the  effective  damping  a  [see  (IV-3)I  is  given 

by 


a 


e 


I  2 

47rL  (jj 


3p 


(1-7) 


where  L  is  the  half-thickness  of  the  film,  p  is  the  resistivity,  and  lo  ^  =  47rM^y.  For  a 

Permalloy  film  with  p  =2X10  abohm  cm,  a  is  equal  to  the  intrinsic  damping  of  about  10~^ 

>  -5  ^  o 

for  L  =  1.6  X  10  cm.  Thus,  unless  the  film  thickness  is  greater  than  about  2500A,  eddy  cur¬ 
rents  may  be  neglected.  Where  necessary,  a  will  be  assumed  to  include  a  ,  but  we  shall  not 
be  concerned  explicitly  with  this  effect  in  the  theory  to  follow. 


5 


We  remark  here  that  we  shall  find  it  useful  to  introduce  the  concept  of  a  "typical  Permalloy 
film"  (TPF),  the  properties  of  which  are  listed  in  Appendix  A.  This  we  do  mainly  for  the  purpose 
of  numerical  estimates  of  our  results,  but  a  few  of  the  approximations  we  use  would  have  to  be 
altered  in  order  to  apply  our  theory  to  a  film  with  parameters  drastically  different  from  a  TPF. 
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IL  MAGNETOSTATIC  FIELD  OF  A  NONUNIFORMLY 
MAGNETIZED  THIN  FILM 


A  fundamental  source  of  difficulty  in  any  magnetic  problem  involving  fluctuations  of  the 

magnetization  M(r)  on  a  scale  comparable  to  a  sample  dimension  is  the  magnetostatic  (dipole- 

dipole)  interaction.  This  long  range  force  leads  to  a  field  H(r)  which  depends  not  only  on 

M(r)  but  also  on  M  throughout  the  sample.^  In  conjunction  with  other  forces  such  as  exchange 

and  anisotropy,  it  usually  renders  static  calculations  very  difficult  (only  a  few  exactly  soluble 
24 

cases  are  known  )  and  dynamic  calculations  intractable.  Thin-film  geometry,  however, 
greatly  simplifies  the  magnetostatic  problem,  enabling  us  to  find  an  approximate  solution  for 
H(r)  in  a  form  which  is  very  convenient  for  both  static  and  dynamic  situations. 

We  calculate  this  field  (also  called  the  stray,  dipolar,  or  demagnetizing  field)  from 
Maxwell's  magnetostatic  equations,  for  a  planar  film  of  thickness  2L  with  an  arbitrary  magnet¬ 
ization  distribution  M(r).  A  rectangular  coordinate  system  is  used  with  film  boundaries  given 

by  X  =  ±L  ,  y  =  ±L  ,  and  z  =  ±L  (z-axis  normal  to  film  plane).  We  eventually  go  to  the  limit 
- 1  ^  ^ 

L  L  —  so  that  edge  effects  may  be  ignored,  i.e.,  the  field  from  poles  on  the  film  rim  is 
neglected.  For  convenience  in  this  section,  we  introduce  a  normalized  position  vector  a  =  r/L 
and  use  the  notations  V  =(  9/9^ ,  9/9t],  9/9?)  and  f'( (7 )  =  9f/9^ .  Maxwell's  magnetostatic 
equations  are  given  by  (I-5a,  b): 


V  .  B((t)  =  0 


(Il-la) 


VX  H((t)  =  0 


(Il-lb) 


where 


B  =  H  +  47rM 

From  (Il-la)  and  (Il-lb)  we  have  the  boundary  conditions 


(II-lc) 


B  •  i  continuous  at  ^  =  ±1 
z 


H  X  i  continuous  at  ^  =  ±1 
z 


(II-2a) 

(II-2b) 


where  i  is  a  unit  vector  in  the  z-direction.  Also,  we  require 
z 


lim  H(a)  =  0 

kl— 


(II-3) 


We  write  M  in  the  form 


M((t) 


m_^(£)  e" 

K 


-1  <  £  <  1 


=  0 


kl  >  1 


(II-4) 


t  In  this  section  we  omit  the  subscript  on  H  for  brevity. 
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where 


(n  n 

IT  ’  IT '  ° 

X  y 

From  (Il-lb),  H  =  -  so  that  the  potential  0(a)  satisfies 

V^0(o‘)  =  47rv  •  M  -  47r  ^  [  i^  •  m_^  (^)  4-  i^  •  m_^(?)]  e^'^'  ^  -  1  <  f:  <  1 


)  (n  ,  n  -  0,  ±1,  ±2,  .  .  .) 
/  X'  y 


-0  |^|>1 

The  problem  now  consists  of  solving  (II-6)  with  the  boundary  conditions 

continuous 


{4n  i^  •  M  —  0  ') 


90 


d(f> 

s-±i'  ^ 


^=±1 

continuous 


^=±1 


and 


lim  0  =  const. 

We  adopt  a  description  of  the  potential  0  of  the  form  (II -4): 

<p(ff)  =  ^  <p_(t)  e^'"'  ^  -I  4  t  <  i 


K 

K  0  for  ^  <  —  1 

Then  from  (II-6)  and  (II-8) 

-K^0_(S)  +  4,'lU)  =  4np^(t) 

K  K  K 


+  0l"(?)  =  0 

K  K 


where  k  =  |  k  |  and 

4  iK  ■ 

K  K  K 

(the  volume  pole  density).  The  transformed  boundary  conditions  are 

47r  i  •  m  (±1)  —  0'  (±1)  =  —0^  (±1) 

^  K  K  7 

<p_^{±i)  -  4>^{±\) 

K  K 

lim  4>^{S)  -  const. 

g— ±oo  K 


(II-5) 


(II-6) 

(II-7a) 

(II-7b) 

(II-7C) 


(II-8) 

(II-9a) 

(II-9b) 

(II-IO) 

(Il-lla) 

(Il-llb) 

(II-llc) 
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For  K  =  0,  the  solutions  are  easily  found  to  be 


(p  =  47r  i 
o  z 


<P^  -  4tt 


r- 

Jo  o 

.  r'  -* 

1  \  m 

z  Jn  o 


du 


(t)  dr 


<t>  =  47r  i 

o  z 


■I 


-1 


m^lf)  df 


(II-12b) 


(II-12b) 


(II-12C) 


to  within  an  arbitrary  additive  constant.  For  k  ^  0,  the  solution  to  (II-9a)  is 

and  from  (II-9b)  and  (II-llc) 


<b  -A  e  +B  e  +  — 

K 

K  K  K 


sinh[K{r  — 

■"0  7 


{II-13a) 


{II-13b) 


K  K 


Using  the  boundary  conditions  (II- 11a)  and  (II- 11b),  the  constants  A_^,  and  C_^  are  found  to  be 

K  K  K 


A 


B  = 


r5h 

1 

II 

C 

.  Ztt  -kt*  -♦  ... 

e  dr  +  —  e  1  •  m  (1) 

K  Z 

K 

1 

II 

C 

K^r  jw  Ztt  ~ k  t*  ^  j  » 

e  dr - e  1  m(—  1) 

K  z 

K 

Ztt 

^-1  K 

Ik  k  J 

so  that 


and 


K 


=  -  ? 


r  P_^(^')  e  ^  ^  ^  ^  ^  di^  —  I 
J-1  7  I 


i  •  m  (i^)  e 
z 


k-r  I 


1 

v  =  -i 


C  p  (r)  dr  -[T  • 

—  .  .  iKl^l 

m  (t')  e  I 

1 

‘J-1  K  1  ^ 

K 

i^  =  -l 

(II-14a) 


(II-14b) 


A  general  answer  to  the  question  ’’Given  the  magnetization  distribution  in  a  thin  film,  what  is 
the  magnetostatic  field?”  is  provided  by  (II-14a,  b),  since 


H  =  h_^(r)  ^  -  v0 


(II-15a) 


with 


h_^(r)  =  — iK0_^(r)  —  i  0^(^) 


(II-15b) 
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inside  the  film,  and  the  analogous  expression  outside  the  film.  However,  in  the  eases  we  shall 
be  dealing  with,  M(a)  is  not  given  but  must  be  determined  statieally  (by  either  minimizing  the 
total  energy  or  requiring  that  the  torque  =  M  X  vanish  everywhere)  or  dynamieally 

(from  the  equation  of  motion  (I-Z)  for  M,  namely  9M/9t  =  potential  (11- 14a) 

will  then  lead  to  an  integro-differential  equation  for  m_^(^)  whieh  might  be  solved  by  approximate 

K 

methods.  However,  a  mueh  simpler  proeedure  —  one  whieh  is  readily  applieable  to  a  large  elass 
of  problems  and  through  whieh  eonsiderable  physieal  insight  ean  be  obtained  —  is  to  assume  that 
m_^  is  independent  of  f .  This  assumption,  whieh  we  eall  the  thin-film  approximation  (TFA),  is 

K 

good  providing  the  Fourier  eomponents  m_^(f )  are  slowly  varying  aeross  the  film  thiekness.  It 

K 

will  be  shown  later  that  for  thin  enough  films,  eomponents  for  whieh  the  TFA  is  not  appropriate 
have  very  high  exehange  energies  and  are,  therefore,  only  weakly  exeited. 

Using  the  TFA  in  (II-14a,  b),  we  find 


<PJ^) 

K 


4tt 

K 


—  (1  “  e  ^  eosh  k^)  + 

K 


i  e  ^  sinh  kK 
z 


(II-l6a) 


—  e  sinh  K  [  —  i  —  ±  1  I 

K  \  K  Z/ 


(II-l6b) 


The  Fourier  eomponents  of  the  magnetostatie  field  are  then  found  from  (II-12a-e),  (H-15a,  b), 
and  (II-l6a,  b)  to  be 


h  =  — 47r  i  i  •  m 
o  z  z  o 


(II-17a) 


and 


h  (^:)=— 47rl—  I—  (1— e  ^  eosh  k^:)  +  i  i  e  ^  sinh 


(k^O) 


inside  the  film,  and 


h"^  =  0 
o 


+  i  e 
z 


(4 


sinh  kK  4  i  eosh 
z 


= —4tt  e^^^sinh/c 


[k{k  *  “z)+  k  'z)j 


(II-17b) 

(II-17c) 

(II-17d) 


above  and  below  the  film.  Note  that  (ll-17a)  and  (ll-17c)  are  also  valid  exactly,  i.e.,  for 
o  o 

The  field  inside  the  film  is  of  primary  interest,  sinee  only  this  field  eontributes  to  the 
torque  on  M.  The  TFA  magnetostatie  torque  may  be  expanded  in  a  Fourier  series 


t(ct)  =  M(^,n)  X  h(ct)  =  Z  ^ 


(ll-18a) 
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where  the  components  of  T  are  given  by 


t(^)  =  m  xh  +  y  m  xh(e:) 

“K  K 

K^O 


t(f)=m  xh(^)  +  m  Xh  +  ).  m  Xh  (?) 

K  K  K  K  -K' 


(II-18b) 


(II-18c) 


r  '  If  * 


'  0 


In  keeping  with  the  TP'A,  it  is  the  ? -average  of  the  torque  components  which  is  of  physical 
interest.  Since  the  h_^(?)  are  the  only  ? -dependent  factors,  we  may  replace  h  ^(?)  by 


h,  =  4  C  h_^(?)  d?  . 

K  K 

Then  from  (II-17b)  we  have  our  final  result  for  the  magnetostatic  field  components 


h  =  — 47r 


x('')  +  i2i2X('')|  •  m_ 


where 


x(k)  =  k  ^  e  ^  sinh / 


x(k)  =  1  - 


(11-19) 

(II-20a) 

(II-20b) 

(II-20c) 


Note  that  since  x(9)  =  1  and  x(0)  =  9,  (Il-20a)  holds  for  k  =  0  also. 

We  next  examine  h_^  in  the  limits  k  »  1  (wavelength  small  compared  with  film  thickness) 

K 

and  K  «  1  (wavelength  long  compared  with  film  thickness).  For  the  ^irst  limit, 


so  that 


x(k)  =  1-  ^  +  .. 


(II-21a) 


(Il-21b) 


and 


h_ 

K 


K 


{k  »  1) 


(II-2lc) 


The  first  term  of  (11-21  c)  is  the  usual  infinite  medium  magnetostatic  field  from  volume  poles  and 
follows  immediately  from  (Il-la-c)  with  periodic  boundary  conditions.  The  second  term  is  a 
small  correction  due  to  the  fact  that  surface  poles  do  not  quite  cancel  themselves  out  for  a  film 
of  finite  thickness. 

For  the  second  limit. 


x(k)  =  1  -  k  +  I  +  .  .  .  (ll-22a) 

x(k)  =  k  -  j  +  .  .  .  (Il-22b) 
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so  that 


K 


^  —4tt 


+ 


1  1 
z  z 


{k  «  1) 


The  first  term  of  (II-22c),  written  with  non-normalized  wave  vector  k  =  ~k/L.,  is 


(II-22c) 


h_^  --4nL  ^  •  m_^  (11-23) 

k  k 

This  is  just  the  planar  magnetostatic  field  in  the  long  wavelength  limit,  calculated  by  some 
26 

authors.  It  is  the  infinite  medium  magnetostatic  field  from  volume  poles  attenuated  by  a 
factor  kL  «  1;  the  reason  for  this  attenuation  is  simply  the  absence  of  volume  poles  from  the 
regions  above  and  below  the  film.  The  second  term  in  (II-22c)  is  the  normal  demagnatizing 
field  from  surface  poles. 

To  summarize,  the  magnetostatic  field  components  are  given  in  the  TFA  by 


h  = 

—  47r  i  i  • 

m 

o 

z  z 

o 

h_  = 

x(kL)  +  T  T  x(kL)l 

■ 

k 

1  k 

J 

k 

(II-24a) 

(II-24b) 


where 


H(?)  =  h  +  )'  h  (II-24C) 

o  L.  ^ 

k^O 

and  X  and  x  are  given  by  (II-20b)  and  (II-20c).  Equation  (II-24a)  is  the  uniform  demagnetizing 
field  of  an  oblate  spheriod  in  the  zero  thickness  limit.  The  first  term  in  (II-24b)  is  from 
volume  poles,  with  an  attenuation  factor  x  due  to  the  presence  of  finite  boundaries.  The  second 
term  in  (II-24b)  is  from  the  surface  poles,  with  an  attenuation  factor  x  due  to  finite  wavelength 
and  therefore  some  cancellation  of  these  poles.  We  now  have  the  thin-film  magnetostatic 
field  in  a  simple  form  which  may  be  used  to  calculate  the  static  magnetization  configuration 
(Sec.  Ill)  and  interactions  in  large-angle  magnetization  rotations  (Sec.  V),  and  as  a  by-product 
we  can  also  obtain  the  planar  mode  spectrum  (including  exchange  fields)  for  all  wavelengths 
(Sec.  V). 
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IIL  STATIC  EQUILIBRIUM:  MAGNETIZATION  RIPPLE 


INTRODUCTION:  EXPERIMENTAL  BACKGROUND 


To  find  the  transient  response  to  an  external  field  of  the  magnetization  IVl(r,  t)  of  a  ferromagnet, 

the  initial  state  M(r,  O)  must  be  known,  which  in  our  formulation  means  specifying  the  initial  val- 

ues  of  the  Fourier  components  m^(t)  for  all  k  .  The  most  important  sources  of  disturbances 

k 

m^(0)  (k  0)  are  thermal  agitation  and  inhomogeneities.  Thermally  excited  spin  waves,  which 
k 

are  present  in  any  ferromagnet  at  a  finite  temperature,  and  which  result  in  the  well-known  de¬ 
crease  of  saturation  magnetization  with  increasing  temperature,  will  be  important  dynami¬ 
cally  if  they  become  unstable.  However,  it  will  be  shown  in  Sec.  V  that  while  such  instabilities 
Z7 

do  occur,  their  growth  times  are  always  longer  than  the  time  for  coherent  switching,  so  that 
they  have  no  significant  first  order  effects. 

Inhomogencities  are  present  to  some  extent  in  any  ferromagnetic  crystal.  However,  in  a 
polycrystallinc  sample,  in  which  we  are  interested  in  this  work,  particularly  large  dispersive 
effects  are  to  be  expected,  both  from  crystalline  anisotropy  forces  and  from  strains  between 
crystallites.  In  this  section  we  calculate  the  equilibrium  dispersion-indueed  components  m  ^(0) 

[henceforth  the  argument  (0)  will  be  dropped]  for  a  polyerystalline  thin  film  (<2000 A),  repre¬ 
sented  by  a  two-dimensional  model.  (The  model,  however,  is  readily  capable  of  extension  to 
three  dimensions  to  treat  bulk  polyerystals  and  thick  films.) 

The  solutions  we  will  find  are  not  only  of  importance  in  specifying  initial  values  in  the 
dynamic  problem,  but  are  also  of  interest  in  their  own  right.  Before  proceeding  with  the  model, 

it  will  be  useful  to  review  the  pertinent  experimental  evidence  concerning  the  equilibrium  mag- 

28 

netization  configuration  in  thin  films.  In  i960  Fuller  and  Hale,  using  a  defoeused  electron 

microscope,  observed  the  Lorentz  deflection  of  electrons  by  the  magnetic  induction  field 

B  47rM(r),  and  discovered  a  wave-like  magnetic  fine  structure  which  they  called  magnetiza- 

29  30 

tion  ripple.  This  discovery  has  been  confirmed  by  subsequent  investigators,  *  and  an  ex¬ 
ample  is  shown  in  Fig.  2. 

Fuller  and  Hale  described  two  basic  ripple  components:  longitudinal  (LMR),  in  which 
M  =  IVl(x)  (k  II  r^),  and  lateral  or  transverse  (TMR),  in  which  M  =  M(y)  (k  ||  i^),  where  the 
mean  magnetization  m^  lies  along  the  x-axis  of  a  Cartesian  coordinate  system  with  the  z-axis 
normal  to  the  film  plane.  These  two  ripple  components  are  shown  schematically  in  Fig.  2.  Since 
the  volume  pole  density  is  very  small  for  LMR  (V  •  M  ^  0)  but  not  for  TMR  (V  •  M  ^  0),  and 
since  the  two  are  equivalent  with  respect  to  exchange,  anisotropy,  and  Zeeman  forces,  the  main 
contribution  to  the  magnetization  ripple  should  be  LMR.  Fuller  and  Hale  thus  interpreted  the 

fine-strueture  lines  they  observed  as  the  loci  of  constant  M,  orthogonal  to  m  .  They  found  a 

-4  -2  ^ 

mean  wavelength  of  --2  X  10  cm  and  a  mean  amplitude  of  --10  radian. 

29 

Fuchs  has  suggested  that  the  origin  of  ripple  is  to  be  found  in  crystalline  anisotropy 
forces,  which  vary  randomly  in  direction  from  crystallite  to  crystallite  (crystallite  size 

^  cm;  crystalline  anisotropy  energy  0^  erg/cm^).  The  magnetization  does  not  follow  these 
local  wanderings  of  the  direction  of  minimum  anisotropy  energy  but,  because  of  exchange  coupling, 
which  tends  to  straighten  the  path  of  M,  follows  the  mean  easy  axis  averaged  over  a  number  of 
crystallites.  For  TMR,  this  number  is  greatly  increased  by  magnetostatic  coupling,  so  that  the 
amplitude  of  TMR  is  much  less  than  that  of  LMR.  The  mean  magnetization  direction  is  deter¬ 
mined  by  the  uniaxial  M -induced  anisotropy  and  the  external  field. 
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In  addition  to  crystalline  anisotropy,  there  will  also  be,  as  just  mentioned,  random  local 
magnetostrictive  forces  due  to  inhomogeneous  strains  between  crystallites  (or  clumps  of  crys¬ 
tallites).  Stresses  (for  example,  surface  tension)  between  adjacent  crystallites  on  the  substrate 
to  which  the  film  is  bonded  result  in  planar  strains;  z-directed  stresses,  on  the  other  hand,  may 
be  for  the  most  part  relieved,  since  the  free  upper  surface  of  the  film  imposes  no  comparable 
constraint.  Thus  we  expect  the  local  magnetostrictive  anisotropy  axes  to  lie  nearly  in  the  film 

plane  and  to  be  randomly  oriented  in  this  plane.  Experimental  evidence  for  such  an  anisotropy 

3 1 

is  provided  by  oblique-incidence  Permalloy  films,  which  are  deposited  by  a  vapor  beam  at  a 
fixed  angle  0^  to  the  substrate  normal.  (In  normal  films  0^  =  0.)  In  such  films,  crystallites 
form  chains  perpendicular  to  the  incident  beam  by  means  of  a  self-shadowing  mechanism,  so 
that  strain  axes  are  not  random  but  tend  to  lie  along  these  chains.  As  a  result  of  both  aniso¬ 
tropic  strain  and  shape  anisotropy  of  the  chains,  these  films  possess  a  large  over-all  uniaxial 
anisotropy  ^  lO^erg/cm^,  depending  on  0^^  and  film  composition.  Shape  effects  may  be  sepa¬ 
rated  from  strain  effects  by  considering  the  compositional  dependence  of  K  ;  from  this,  one 

4^/3 

may  estimate  random  anisotropy  energy  in  normal  films  of  ~5  X  10  erg/cm  for  a  TPF. 

The  relative  roles  played  by  crystalline  and  strain  anisotropies  in  inducing  magnetization 
ripple  have  not  been  firmly  established  experimentally.  In  what  follows  we  shall  neglect  crys¬ 
talline  anisotropy  since  strain  effects  are  probably  greater  and  strain  anisotropy  is  easier  to 
handle  analytically.  However,  consideration  of  crystalline  anisotropy  would  not  alter  the  results 
in  any  significant  qualitative  way.  We  assume  a  random  local  magnetostrictive  anisotropy  which 
is  planar,  uniaxial,  and  constant  across  the  film  thickness  and,  using  a  simple  model  for  a  poly¬ 
crystalline  film,  proceed  to  calculate  the  ripple  spectrum  in  static  equilibrium. 

B.  RANDOM  ANISOTROPY  MODEL 

Our  model  for  a  film  is  a  planar  array  of  N  right  cylindrical  cells  of  height  2L  (=  film 
thickness)  and  cross-sectional  areas  a  (m  =  1,  2,  ,  .  .  N).  The  a  are  assumed  to  be  random 
variables  with  ensemble  average  a,  and  the  cross-sectional  shapes  are  also  assumed  to  vary 
randomly.  By  this  we  mean  precisely  that  the  intersections  of  any  line  in  the  film  plane  with 
the  cell  boundaries  occur  at  random  (uncorrelated)  intervals  along  the  line  (see  Fig.  3). 


U-?4-5S>0| 


Fig.  3.  Rondom  ceil  model  showing 
construction  used  to  find  outocorre- 
lotion  function  c  (  r  ). 


o 


14 


Furthermore,  the  mean  number  of  such  intervals  per  unit  length,  r^  ^ ,  is  assumed  to  be  the  same 
for  all  lines  in  the  plane;  i.e.,  cell  shapes  are  on  the  average  isotropic  and  cell  sizes  and  shapes 
are  on  the  average  homogeneous.  The  film  plane  is  of  area  S  =  Na,  and  we  shall  eventually  pass 
to  the  limit  N  «  (a  fixed). 

Each  cell  has  a  magnetic  anisotropy  which  is  the  sum  of  two  terms;  a  uniaxial  M-induced 

anisotropy,  common  to  all  cells,  with  energy  density  and  easy  axis  the  x-axis,  and  a  uniaxial 

strain-induced  anisotropy  with  energy  density  and  axis  in  the  film  plane  at  an  angle  to 

the  x-axis,  where  may  take  on  any  value  between  —  tt/Z  and  tt/z  with  equal  probability.  We 

further  assume  zero  correlation  among  the  quantities  a^,  and  and  that  the  variations 

in  and  are  uncorrelated  from  one  cell  to  the  next.  This  addition  of  uniform  and  random 

uniaxial  anisotropies  has  been  used  before  to  explain  some  film  properties, but  a  model  was 

lacking  and  exchange  and  magnetostatic  interactions  were  neglected;  in  the  present  work  inter- 
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actions  play  a  vital  role.  Rother  has  analyzed  a  model  similar  to  ours;  however,  in  his  model, 

=  0  and  the  cells  have  equal  square  cross  sections.  He  has  included  the  effects  of  exchange 
interactions,  and  magnetostatic  interactions  as  a  perturbation,  but  his  results  are  quite  different 
from  ours  for  reasons  which  will  be  made  clear  at  the  end  of  this  section. 

The  anisotropy  energy  density  at  any  point  in  the  film  is  given  by  an  obvious  generalization 
of  (1-6) 

E  (?)  = -M  {K  [M(r  )  •  T  ]^  +  K.(r)  [M(r  )  •  tt(r  )]^}  (III-l) 

a  o  o  ‘  x'  1  ‘  ^  ^ 

where  a(r  )  is  a  unit  vector  at  angle  a  to  the  x-axis;  for  r  in  the  m^^  cell,  K^(r  )  -  and 

air  )  =  The  anisotropy  field  is  then 

H  =  -  V  E  =  2M  [K  (T  -  mT  •  MM  T  •  M  +  K.i?  -  M?  •  MM  a  ■  Ml 
a  '  — ^  a  o'o'x  X  ox  1'  o  ^ 

so  that  the  anisotropy  torque  is  given  by 

f  =  M  X  H  =  2M  "^M  X  (K  T  T  +  K.aa)  ■  M  .  (III-2) 

a  a  o  o  X  X  1 

Since  we  shall  be  interested  only  in  the  ease  of  external  fields  in  the  film  plane,  M  must  lie  in 
this  plane  in  static  equilibriumt  and  is,  therefore,  determined  by  the  angle  (p{r  )  between  M  and 
the  x-axis.  Equation  (III-Z)  then  becomes 

T  =  i  [— K  sin  Z(^  +  K.  sin  Z(cv  —  )]  .  (III-Z-1) 

a  z  ‘  o  1 

We  next  separate  the  spatially  varying  part  of  (p,  writing 

(p{r)  =  +  d(p{r)  (III-3a) 

where 

<Pq  =  ~  ^  ^  (ni-3b) 


and  from  (I-la) 

|6<^  I  «  1  .  (III-3c) 


t  This  moy  not  be  true  if  there  is  a  strong  eosy  oxis  normol  to  the  film  plone;  but  in  thot  cose,  deviofions  of  M 
from  its  meon  direction  ore  lorge,^^  ond  thus  beyond  the  scope  of  this  lineorized  treotment. 
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(III-4) 


2 

Inserting  (III-3a)  in  (III-2-1)  and  dropping  terms  of  order  {6(p)  and  higher,  we  find 

T  =T  (  — K  sin  2<p  +  P(r  )  —  26</?(r  )  [K  cos  2</7  +  Q(r  )]) 

a  z  o  o  ^  ^o  ^o  ^"^ 

where 

P(r  )  =  K^(r  )  sin  2  [a(r  )  -  <p (III-5a) 

Q(r  )  =  K^(?)  cos  2  [^(r  )  —  (/?^]  .  (III-5b) 

The  functions  P  and  Q  contain  all  effects  of  the  random  anisotropy  field.  P  may  be  thought  of 
as  a  random  perturbing  force  on  6<p,  and  Q  as  a  random  restoring  force.  The  ensemble  aver¬ 
ages  of  these  functions  vanish,  since 

P  =  sin  2(q'  —  (p^) 


=  K.  sin  2(q'  —  (P  )  (since  K.  and  a  are  uncorrelated) 
1  ^o  Im  m 

A  p  7r/2 

=  K.  —  \  sin  Z(a  —  (p  )  da  -  0 

^  ^  J-tt/Z 


and  similarly  for  Q. 

We  now  expand  dcp,  P,  and  Q  in  Fourier  series 

6<p{r  )  =  ^  (p^  c^^  ’  ^ 

k  7^0  ^ 

P(r)= 

k 

Q(r  )  =  E  Q-* 

7*  k 


(III-6a) 


(III-6b) 


(III-6c) 


where  the  wave  vector  k  takes  on  values  imposed  by  periodic  boundary  conditions  at  the  film 
edge  X  =  y  =  ^^y* 


k 

with  the  condition 


(n  ,  n 
X*  y 


0,  ±1,  .  .  .  ) 


(III-7) 


ki  «  1 


(i  =  lattice  constant)  insuring  the  validity  of  our  continuum  treatment  of  the  magnetization.  The 
components  of  (III-4)  then  become 


(III-8a) 


(III-8b) 
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These  expressions  will  be  eombined  with  the  remaining  torques  aeting  on  M  (exehange,  mag- 
netostatie,  and  Zeeman),  to  determine  the  ripple  eomponents  (p .  In  partieular,  we  will  be 

- 2  ^ 

interested  in  the  "power”  spectrum  |<^_^|  ,  and  to  find  this  ensemble  average  will  require  the 
- 2  ^ 

power  speetrum  |  f _^  |  of  a  random  funetion  F  of  the  form 


F(r  )  =  F 


for 


.  th  ,, 
r  in  m  ecll 


(III-9a) 


F(r  )  =  0 

where  F  may  be  P  or  Q,  and 


(III-9b) 


F(r  )  =  ^  e 


ik  •  r 


(k  in  film  plane) 


(III-9e) 


To  aceomplish  this  we  first  find  the  autoeorrelation  function 


Co(r  ) 


F{r’  +  r  )  F(r')  d  r' 


(III-IO) 
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by  a  generalization  of  a  method  due  to  Kcnrick.  Consider  a  straight  line  in  the  film  plane  de¬ 
fined  by  the  veetor  r,  with  origin  at  an  arbitrary  point  r'  (see  Fig.  3).  We  first  determine  the 
probability  PqCi',  r')  that  no  eell  boundary  interseets  this  line  between  r'  and  r'  +  r  .  From 
homogeneity  and  isotropy,  respectively,  ^  probability  that 

there  are  no  eell  boundaries  between  0  and  r  +  Ar  is  equal  to  the  probability  that  there  are  none 
between  0  and  r,  multiplied  by  the  probability  that  there  are  none  in  the  interval  Ar.  This 


latter  probability,  in  the  limit  Ar 
along  any  line).  Thus 


0. 


-1 , 


is  1  —  r  Ar  (where  r 
o  o 


-1 


is  the  density  of  cell  boundaries 


lim  P  (r)  (1 
Ar-O  ° 


r  Ar)  =  lim  P  (r  +  Ar) 

^  A  n  O 

Ar-^0 


P  (r)  +  ArP' (r)  +  . 
o  o 


or 

p;(r)  =  -rj^p^(r) 


which  has  the  solution 


P 

o 


(r)  =  e 


-r/ r 


o 


(III-ll) 


[since  Pq(0)  = 

Since  F(r  )  is  a  random  funetion  with  zero  mean,  the  integrand  of  (III-IO)  vanishes  unless 
r’  +  r  and  r’  are  in  the  same  cell,  in  which  case  it  is  equal  to  F^.  Thus  we  have  the  result 


^  — 7  -r/r 

e^(r  )  =  e^(r)  =  F  c  ^ 


(III-12) 


This  expression  is  somewhat  in  error  for  two  reasons.  First,  it  ignores  the  effects  of  finite 
boundaries,  which  would  result  in  a  correction  of  order  N”^  and  may  therefore  be  neglected. 
Second,  it  does  not  consider  the  possibility  that  because  of  a  concave  eell  boundary,  both  r 
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and  r’  may  lie  in  the  same  cell  although  the  straight  line  connecting  them  passes  through  one 
or  more  other  cells  (as  in  cell  m^.  Fig.  3).  However,  unless  the  cells  are  serpentine,  which 
is  not  at  all  plausible  physically  since  crystallites  grow  outward  in  all  directions  from  a  nucleus, 
contributions  from  such  double  crossings  (and  triple,  etc.)  should  be  very  small  and  can  safely 
be  ignored. 

The  power  spectrum  is  just  the  Fourier  transform  of  the  autocorrelation  function: 


1 2  If  -ik.  r  ,2-> 

.1  =  s  j  'I  ^ 


S 


-r/r 

r  e  ^  rdr  C 


Ett 


-ikr  cos  0 
e  do 


exp 


/  2  2 
V  X  +  y 


—  i(k  X  +  k  y) 
X  y*^ 


dxdy 


(III-13) 


where  S  is  the  region  |x|  >  L  ,  |y|  >  L  .  The  first  term  in  the  bracket  in  (III-13)  is 

X  y 

-r/ 


Ztt  r  e  °  J  (kr)  rdr  =  27rr  ^(1  +  k^r 
1  o  o  o 

o 


The  second  term  vanishes  as  L  ,  L 

X  y 


exp 


/  2  ,  2 

—  — ^ —  i(k  X  +  k  y) 
r  X  v*^ 


dxdy 


^  ^ \  exp 
S 


/  2  ~  2 
NX  +  y 


dxdy 


r 


-r/r 
e  ^  rdr  C 
*^0 


Zti 


d0  =  27rr  (L  +  r  )  e 
o  s  o 


-L  / 


s  o 


which  goes  to  zero  as  L  where  L  is  the  smaller  of  L  and  L  .  Thus  we  have  the  main 

®  s  s  X  y 

result  of  this  section 


^  27rr^^(l  +  .  (III-14) 

It  may  be  helpful  to  know  the  mean  cell  boundary  spacing  r^  in  terms  of  the  cell  areas  a^; 
this  is  most  easily  found  by  direct  evaluation  of  (f^l  ^  comparison  with  (III-14).  The  k  =  0 
component  of  F  is 

=  I 

N 

^  "S  S 

m  =  l 


so  that 


18 


giving 


N  N 


ifi^  =  ^y  vff 

'  o  ^  t-i  m 

^  m=l  n=l 


_  3.  3. 

m  n  m  n 


F^  V  2 

c;,2  L 

^  m  =  l 


(since  F  F  =  F  6  ) 

m  n  m,  n 


1:^2  2 
F  3 

s  - 

3 


(since  N  =  S/3) 


a" 


\  Zna.  J 


Fin3lly,  from  (111-53,  b). 


.2  .2 


1  ..2 


P  =  sin  2  [Q!(r)  -  (p ^  =  Q 


(III-15) 


(III-16) 


C.  MAGNETIZATION  RIPPLE  SPECTRUM 

The  rem3ining  contributions  to  the  torque  3re  now  collected.  We  consider  the  C3se  of  3 
ste3dy,  uniform,  extern3l  magnetic  field  H  in  the  film  pl3ne  at  an  angle  (3  to  the  x-axis  (see 
Fig.  4).  The  Zeeman  torque,  in  the  linear  approximation,  is  then 

f ,  =  M  X  H 
h 


=  i  M  H[sin(/?  — )  —  (5</?  cos{S  —  (p  )1 
zo  ^^o  ^o^ 

with  Fourier  components 

( t  ),  =  i  M  H  sin  {(3  —  <p  ) 

'oh  z  o  ^  ”0 


(III-17a) 


(III-17b) 


^  cosO  -  <p^)  (p_^ 


(III-17C) 


Fig.  4.  Coordinate  system  for  static  case. 


|3-?4-55ai| 
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The  exchange  field  is  given  by  (1-4) 


H  = - M(r)^^^(— i  sine/?  +  i  cos(p  )V  6(p 

e,-2  M  x  y  ^ 

Mo 


(III-18a) 


where  A  is  the  exchange  constant.  This  field  exerts  an  exchange  torque 

T  =  M  X  H  =  T  2AV^6w 
e  e  z  ^ 

with  Fourier  components 
=  0 

o  e 


i  2Ak  (p 

k  'e  ^  k 


(III-18b) 

(III-lSc) 

(III-18d) 


Finally  we  have  the  magnetostatic  field  with  Fourier  components  in  the  TFA  given  by 
(II-24a,  b),  which  become 


(h  )  =0 

o  m 


(III-19a) 


k  k  •  m 


-47rx(kL) - ^ 

k'm  k 


(III-l9b) 


where 


X  (k  /  =  1  —  K  e  sinh  k  =  k  —  k  +  .  .  . 


(in-20) 


The  linearized  torque  components  are  then 


( t  )  =0 

o  m 


-*  2r^  2 

=  “  i  47rM  x'(kL)  sin  —  <p  )  <p^ 
k'm  k 


(III-21a) 

(III-21b) 


where  k  lies  in  the  film  plane  at  an  angle  ^  to  the  x-axis  (see  Fig.  4). 

The  equilibrium  value  of  6<p  is  determined  by  the  condition  that  the  effective  torque  vanishes, 
or 

ik  •  r 


T  .  J  r  )  =  T  +  T,  +  T  +  T 
eff  ahem 


m 


=  0 


(III-22) 


k' eff 


The  Fourier  components  of  (III-22)  vanish  independently,  so  that  from  (III-8a),  (III- 17b), 
(III-18C),  and  (III-21a) 


—  K  sin  2(p  +p— 2  V  q  w  +MH  sin(fl  —  w  )  =  0 

O  ^0^0  LJ  ^  7^  r*  O  O 


(III-23a) 


ky^O 


-k  k 


Fork  0,  from  (III-8b),  (III-17c),  (III-18d),  and  (III-21b) 
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(III-23b) 


—  2K  Qos2(p  (p  +  p  —  2 
®  °  k'  k 


Z 

k'^O 


q  (p  —  [  M  H  cosip  —  (p  )  +  2Ak^ 
k-k'  k'  °  ° 


2r^  2 

+  47rIV1^5r  ^ 


\Ve  solve  (III-23a,  b)  by  an  iterative  procedure,  first  assuming  (p -  0  for  all  k  ^  0  and  solving 
for  (p Equation  (III“23a)  becomes  ^ 


1 

2 


sin  2(/7^  —  h  sinifi  —  (p 


K 


(III-24) 


where  the  reduced  field  h  is  given  by 


h  = 


H 


K 


HM 

_ q 

2K 

o 


(III-25) 


The  term  p^  is  the  magnitude  of  the  (spatial)  average  torque  on  M  due  to  random  anisotropy 
forces,  and  is  only  nonzero  because  for  any  given  ensemble  of  values  of  a  ,  a  ,  and  K. 
there  is  a  small  but  finite  probability  that  this  average  does  not  vanish.  Of  course  the  ensemble 
average  p^  =  0. 

The  mean  square  value  of  the  right  side  of  (III-24),  obtained  with  the  help  of  (III-14)  and 
(III-16)  is 


8  S 


(III-26) 


where  V  is  the  ratio  of  the  rms  value  of  the  random  anisotropy  energy  to  the  uniform  anisotropy 
energy,  or 


r 


2 


(III-27) 


Since  from  (III-26)  the  right  side  of  (III-24)  vanishes  in  the  limit  S  we  may  neglect  it  in  a 

first  approximation  to  (p obtaining 

2  sin  Z(p^  —  h  sin  {(3  —  (p  =  0  (III- 28) 

1 3 

which  is  the  well-known  equilibrium  equation  for  the  magnetization  direction  of  a  uniaxial  film. 
We  shall  be  interested  only  in  three  special  eases  of  this  equation,  for  which  solutions  are 

(P^  =  0  (^  =  0,  h  >  ~1)  (III-29a) 

<p^  =  tt/z  (/3  =  tt/Z,  h  >  1)  (III-Z9b) 

(p^  =  arc  sinh  (/?  =  tt/z,  |h|  <  1)  .  (III-Z9e) 

It  is  easy  to  see,  by  examining  the  average  uniaxial  anisotropy  and  Zeeman  energies 
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(III-30) 


2 

E  =  K  sin  (D  —  M  H  cos  —  (p  )  , 

O  O  ”0  O 

that  the  solutions  (III-29a-c)  minimize  E^,  and  are  therefore  stable. 

The  next  step  in  the  iteration  procedure  is  to  assume  that  there  is  but  one  nonvanishing 

^  0),  and  using  the  value  we  have  just  obtained  for  (p  ,  solve  (III-23b)  for  tp  ,  with 
k  ^  ^ 

<p_^{k^  ^  k,  0)  set  equal  to  zero.  This  just  means  that  to  a  first  approximation  the  ripple  com- 
k' 

ponents  cp are  decoupled.  (Note  that  it  is  only  through  the  random  anisotropy  components 
k 

q  that  (p  and  (p  are  coupled.)  Equation  (III-23b)  becomes 
k-k'  k  k’ 


+  cos  2(p 
K  ^  o 

o 


+  h 


cos  {p  —  (p^) 


Ak 

K 


27rM  ' 
c 

K 


X  sin 


■(4.  -cp^) 


K 


(III-31) 


4  _  ^ 

The  term  q^KJ^  in  (III-31)  may  be  safely  neglected,  for  the  same  reason  that  the  term 
was  negligible  in  (III-24),  provided 


A(/3,  h)  =  cos2(/7^  +  h  cosip  —  (p ^  0  (III- 32) 

with  (p^  determined  by  (III-28).  But  A  is  proportional  to  d(p ^  [see  (III-30)]  so  that  for 

stable  equilibrium 


A  >  0 


(III-33) 


and  (III-32)  is  satisfied.  The  set  of  values  ip^,  h^)  such  that 


A(/?^,  h^)  =  0 


(III-34) 


occurs  at  transitions  from  stable  to  unstable  equilibrium  and  are  thus  threshold  angles  and  fields 
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for  irreversible,  coherent  rotation  of  m^.  We  may  expect  serious  difficulties  in  the  theory 
near  such  thresholds. 

For  A  0,  (III-31)  becomes 


P_^ 

(p  =  o  T? — 7  [1+Xk  +X  L  Y  sin  (4>  —  (Z?  )  1 
2KA^  e  m  ^o' 

k  o 

(ky^O) 

where  we  have  defined  an  "exchange  length" 


and  a  "magnetostatic  length" 
27rM 

X  =  _ ^ 

'^m  K  A 
o 


(III-35) 


(III-36) 


(III-37) 


A(/3,  h),  defined  by  (III-32),  is  for  cases  (a),  (b),  and  (c),  respectively,  of  (III-29) 
A(0,  h)  =  h  +  1  (h  >  -1) 

A>(|,h)  =  h-1  (h>l) 


(III-38a) 

(III-38b) 
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A<(|,h)  =  1  -h^  (|h|  <  1) 


(III-38c) 


The  mean  square  ripple  amplitudes  then  follow  from  (III-14),  (III-16),  (III-27),  and  (III-35): 


k 

(k?^0) 


,,  ,  Z.2.-3/Z  ,  ,-2 

— ?  — g—  ( 1  +  k  )  ( 1  +  g_ ) 


8A 


(III-39) 


where 

g  =\^k^  +  \  L  ^  X  sin^  (4>  —  (^  )  >  0  .  ( III -40) 

k 

Equation  (III-39),  the  ripple  spectrum,  is  the  main  result  of  Part  C.  In  Appendix  B  we  return  to 
{III-23b),  and  carrying  our  iteration  procedure  one  step  further  show  that  the  first  order  coupling 
[meaning  the  (p  are  the  coupling-free  values  of  (III-35)]  introduces  a  negligible  correction  to 
- 2 

•  This  gives  us  some  confidence  in  our  method  of  solution  and  in  the  result  [(III-39)], 
k 

which  we  shall  use  in  Sec.  V.  In  Sec.  III-D,  however,  we  examine  this  result  and  show  that  it 
represents  a  magnetic  fine  structure  with  the  characteristics  of  the  observed  ripple. 

D.  DISCUSSION  OF  RESULT 


First,  it  must  be  pointed  out  that  even  though  the  mean  square  Fourier  components  given 
by  (III-39)  are  proportional  to  S  and  therefore  vanish  in  the  limit  S  ♦  their  density  in  the 
k -plane,  from  (III-7),  is 


S 

iZir)^ 


(III-41) 


Thus  any  observable  effects  of  the  ripple  components  are  independent  of  film  area  S,  provided 
such  effects  do  not  become  too  large  for  long  wavelengths  (of  the  order  of  L  or  L  ).  An  ex- 

^  y 

ample  of  an  observable  effect  is  the  mean  square  magnetization  dispersion 


<5^  =  <[<^(r  )  -  =  <[i5<p(r  )]^> 


(III-42) 


which  we  calculate  in  Sec.  III-E.  From  Parseval's  theorem  [or  from  (III-IO)  and  (III-13)] 


<5^  =  Z  (III-43) 

k^O  ^ 

and  if  the  main  contribution  to  this  sum  is  from  components  with  wavelength  short  compared  to 

L  and  L  ,  we  may  replace  it  by  an  integral  using  the  correspondence 
X  y 

Z  f-.  ^  C  f(k)d^k  (III-44) 

^  k  {Z-Kf  J 


which  follows  from  (III-41).  [f(k)  is  just  f_^  with  k  considered  a  continuous  variable.]  We  see 

2  ^ 
then  that  6  is  independent  of  S,  which  makes  physical  sense. 
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Next  we  consider  the  first  factor  of  (III-39): 


2 


2 


f  2K  cos  2(p  + 

^  o  o 


M  H  cos  {(3  —  (p  )] 


[3^E 


o  ^o 


2,2 


(111-45) 


The  numerator  of  (III-45)  contains  the  perturbing  effect  of  random  anisotropy  forces  (the  factor 
2  coming  from  the  average  over  axis  orientations  cv  ).  The  denominator  may  be  thought  of  as  a 
uniform  magnetic  stiffness  tending  to  align  M  along  11  and  along  the  uniform  easy  axis  (x-axis). 
In  strong  fields  (h  »  1)  the  ripple  amplitudes  go  to  zero  and  the  film  may  be  considered  satu¬ 
rated.  As  the  field  is  lowered  and  a  coherent  threshold  h^  is  approached,  A  0,  the  ripple  am¬ 
plitudes  grow  very  large,  the  condition  (III-3c)  is  violated,  and  we  have  a  threshold  catastrophe. 

The  static  equilibrium  problem  becomes  nonlinear  in  the  Fourier  components  (P^,  and  further 

—  k 

progress  is  very  difficult  without  some  simplified  ansatz  for  M(r  )  such  as  the  band  model  of 

Thomas^^  or  the  labyrinth  model  of  Smith  and  Harte.^^  Such  treatments  are  not  within  the  scope 

of  this  work.  In  the  absence  of  the  uniform  stiffness,  and  in  the  absence  of  the  interaction  term 

,  the  components  cp ^  do  not  of  course  become  infinite:  M(r  )  just  follows  the  random 
k  k 

anisotropy. 

A  crystallite  cutoff  of  the  ripple  spectrum  results  from  the  effects  contained  in  the  factor 

(1  +  (III-46) 

of  (III-39):  for  wavelengths  27r/k  shorter  than  the  cutoff  length  ^Trr^,  ripple  amplitudes  arc 
sharply  diminished.  This  occurs  because  components  of  the  random  anisotropy  with  wavelengths 
shorter  than  the  mean  crystallite  size  are  very  small. 

Interaction  effects,  which  are  contained  in  the  last  factor  of  (III-39) 


k 


(III-47) 


are  of  great  importance  in  determining  the  characteristics  of  the  magnetization  ripple.  The 

term  g  ,  defined  by  (III-40),  is  the  exchange  plus  magnetostatic  interaction  energy  density  for 
k 

unit  |<^>_^|  normalized  to  the  uniform  magnetic  stiffness.  The  effect  of  these  interactions  is 
k 

to  suppress  all  Fourier  components  (p ^  for  which  g_^  »  1,  and  leave  unaffected  those  for  which 

g  «  1.  If  we  define  a  characteristic  closed  curve  C  in  the  k -plane  by  the  relation 
k 

g^  =  1  for  k  on  C'  (III-48) 

k 

ripple  components  with  wave  vectors  on  C  are  attenuated  by  a  factor  of  |  due  to  interactions. 

Since  g_^  monotonically  increases  with  increasing  k,  the  interactions  (roughly  speaking)  suppress 
k 

those  components  with  wave  vectors  outside  C  while  leaving  unchanged  those  with  wave  vectors 
inside  C . 

It  is  helpful  to  transform  C  to  C  in  the  reciprocal  wave  vector  plane,  defined  by  27rk/k^. 
Thus,  points  on  C  are  described  by  a  vector  in  the  k -direction  whose  length  is  the  wavelength 
27r/k.  Figure  5  shows  the  characteristic  curve  C  for  a  TPF  (unless  otherwise  noted,  a  TPF  is 
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Fig.  5.  Characteristic  curve  C 
showing  effect  of  interoctions 
on  ripple  components  with  wove 
vectors  k  (for  o  TPF). 


-5  -  2 

in  zero  external  field),  for  whieh  \  =  2.5  x  10  em  and  \  =  1.25  x  10  em.  In  aeeord  with 

e  m 

our  rough  interpretation  given  above,  all  Fourier  components  of  the  ripple  with  wavelengths  in'- 
side  C  are  suppressed  by  exchange  and  magnetostatic  interactions,  while  components  with  wave¬ 
lengths  outside  C  are  unattenuated.  It  is.  easily  seen  that  the  shortest  wavelength  on  C  is  for 

^  =  <p  ,  or  wave  vector  in  the  direction  of  mean  M.  The  wavelength  of  this  LMR  is  27rX  (=  1.6  X 
_  ^  o  e 

10  em  for  a  TPF),  and  it  is  apparent  from  the  sharpness  of  C  near  the  origin  that  this  ripple 
is  well  oriented,  in  aeeord  with  experimental  observation. 

In  Fig.  6  we  show  this  region  in  greater  detail,  where  now  we  plot  angle  of  wave  vector  to 
mean  M  vs  wavelength  27r/k  on  a  logarithmic  scale.  Equipartition  of  exchange  and  magneto¬ 
static  energies  occurs  along  the  dashed  line,  whieh  satisfies  the  equation 

\^k^  =  \  L"^x(kL)  sin^  (<J>  -  (p  )  .  (III-49) 

e  m  ^  o 

We  see  from  this  that  over  almost  the  whole  reciprocal  wave  vector  plane  inside  C  the  inter¬ 
action  is  predominantly  magnetostatic.  For  wavelengths  much  greater  than  the  film  thickness, 
or 

kL«l  ,  (III-50) 

(III-48)  becomes 

\^k^  +  \  ksin^(4>-</)  )  =  1  .  (III-51) 

e  m  o 

The  condition  (III-50)  is  fairly  well  satisfied  on  C  even  for  the  shortest  wavelength  for 

whieh  kL  =  l/x^(=  0.2  for  a  TPF).  In  the  magnetostatic  region  we  may  neglect  the  first  term 
in  (III-51),  so  that  C  is  determined  by 

(X  k)"^  =  sin^(^“</>  )  .  (III-52) 

m  ^o 

From  this  we  find  that  the  minimum  wavelength  of  TMR  (k  -L  mean  M)  is  ^ttX^  (=  0.8mm  for  a 
TPF). 
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2  TT /k  ( cm) 


Fig,  6,  Detoil  of  Fig.  5  near  origin  showing  exchange-coupled 
ond  mognetostotic-coupled  regions  (far  a  TPF). 


The  characteristic  curve  of  Fig.  5  may  also  be  interpreted  as  the  extent  of  the  interactions. 
M  at  the  origin  is  strongly  coupled  to  M  everywhere  inside  C,  where  Z7rk/k  has  the  meaning 
of  a  position  vector. 

We  have  seen  that  short  wavelength  ripple  components  have  well-oriented  wave  vectors  in 
the  direction  of  m  and  include  all  wavelengths  longer  than  about  ZttX  .  However,  electron 

o  _4  ^ 

micrographs  show  a  fairly  well-defined  wavelength  of  1  to  Z  X  10  cm  (see  Fig.  Z);  one  wonders 

why  this  pattern  is  not  completely  obscured  by  the  longer-wavelength  components.  The  answer 

is  that  what  are  observed  are  spatial  variations  in  the  intensity  of  the  electron  deflection  pattern 

caused  by  spatial  variations  in  the  l^rentz  force.  The  more  rapid  the  variation  in  M,  the 

Z8 

greater  the  electron  contrast.  As  shown  by  Fuller  and  Hale,  the  magnitude  of  the  change  in 

intensity  of  the  electron  beam  due  to  the  k^^  component  of  LMR  is  proportional  to  /|k(/?_^[^  . 
From  (III-39),  ^ 


\k(f>_ 


r  /“  X  2,  2, -3/4  ,  ,,  ,2,  2.-1 

—  N  IT  (l  +  r„k)  k(l+Xk) 


2A 


nTS 


(111-53) 


SO  that  for 


kr  «  1 
o 


(111-54) 


we  find 


\k(p  1^  oc  k(l  +  .  (III-55) 

n/  k 

From  (III-55)  we  see  that  the  intensity  spectrum  has  a  peak  at  k  =  Ordinarily  X^  is  much 

greater  than  the  crystallite  size,  and  (III-54)  is  satisfied.  (For  a  TPF,  r  /\  =  0.05.)  The 

-4  o  e 

wavelength  at  the  peak  is  ^ttX^  (=  1.6  X  10  cm  for  a  TPF),  which  corresponds  quite  closely  to 
the  actual  wavelength  observed. 


Z6 


With  an  external  field  applied  along  an  easy  or  a  hard  axis  of  the  uniform  anisotropy,  cases 
(a)  and  (b)  of  (III-38)  give 

ZttX^  =  [K^(h  ±  for  h  >  =f1  .  (III-56) 


As  the  field  is  increased,  the  peak  wavelength  decreases  until  finally  (III-54)  is  violated.  How¬ 
ever,  this  happens  at  a  field 


(=  800oe  for  a  TPF) 

o  o 


which  is  so  high  that  M  is  nearly  saturated  and  the  ripple  has  disappeared.  Therefore,  it  is 

-5  -4 

only  for  a  film  with  very  large  crystallites  (  >10  to  10  cm)  that  crystallite  cutoff  from  the 
term 


(1 


^  2,  2, -3/4 

+  r  k  ) 
o 


of  (III-53)  need  be  considered.  As  the  field  approaches  a  coherent  threshold  (h  =  t  1 ),  becomes 
infinite  —  another  aspect  of  the  threshold  catastrophe  just  described.  The  actual  conversion  of 
LMR  into  uniform  magnetization  is  prevented  by  nonlinear  interactions  which  limit  the  growth 
of  the  peak  ripple  wavelength.  Similar  considerations  apply  to  case  (c)  of  (III-38),  for  which 

Zt\  =  27rA^/^  [K  (1  -  for  |h|<l  .  (III-57) 

e  ‘  o 


E.  MAGNETIZATION  DISPERSION 

A  vital  aspect  of  magnetization  ripple  remains  to  be  examined,  namely,  its  over-all  mag¬ 
nitude.  Perhaps  the  best  measure  of  this  is  the  magnetization  dispersion  6,  the  rms  angular 
deviation  of  M(r  )  from  m^,  given  by  (III-42).  This  quantity  is  of  great  interest  for  a  number 
of  reasons.  First,  6  must  be  no  smaller  than  the  observed  LMR  if  our  theory  is  to  account  for 
ripple.  Second,  if  our  linearized  treatment  is  to  be  valid  and  (III-3c)  is  to  hold,  we  must  find 

6  «  1  .  (III-58) 


Third,  6  may  be  either  measured  directly  or  compared  with  "anisotropy  dispersion,"  measured 
by  a  variety  of  techniques^  but  always  near  a  threshold  for  coherent  rotation.  Finally,  in 
Sec.  VI  we  will  relate  dynamic  effects  of  the  ripple  to  6. 

From  (III-43)  and  (III-44) 


6^  =  — kdk  r  I<p(k)l^d'i>  (III-59) 

iZn) 

where  |<^(k  )|  is  given  by  the  right  side  of  (III-39).  Replacement  of  the  sum  over  discrete  k 
by  an  integral  is  clearly  valid  here,  since  there  is  no  divergent  contribution  to  the  integral  near 
k  =  0;  for  the  same  reason,  the  fact  that  the  integral  goes  over  k  =  0  while  the  sum  did  not  is 
of  no  importance  whatsoever  in  the  limit  S  — 

From  (III-39),  and  substituting  g_^  from  (III-40),  the  integration  over  is  readily  performed 

.  .  k 

giving 


•  (2  +  2X^k^  +  X^L'^x  )  kdk  (III-60) 
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(III-61) 


_  \ 

We  now  assume  that  the  main  contribution  to  (III-60)  comes  from  k  =  0(X  ).  If 

e 

L  «  \ 

e 

which  is  only  moderately  well  satisfied  for  a  TPF  -  0.2),  we  may  approximate  x(kD  by 

its  first  term,  kL.  [See  Eq.  (III-20) .]  Therefore,  since  for  nearly  all  films  of  practical  interest 

^o  ^e  ^m  (111-^2) 

(for  a  TPF,  =  5  x  10  and  ^^Aj^  =  2X10“  ),  the  integral  in  (III-60)  has  the  approximate 

value  (which  is  also  a  lower  bound  since  x  < 


-1/2 


(1 

e 


(III- 63) 


The  integrand  in  (III-63)  has  its  maximum  value  at  k  =  l/(X^'/^);  this  gives  us  some  confidence  that 
the  approximations  used  in  obtaining  (III-63)  are  good.  A  more  rigorous  approach  —  one  which  is, 
in  fact,  necessary  if  (III-61)  or  (III-62)  are  violated,  is  to  split  the  integral  of  (III-60)  into  seg¬ 
ments  and  in  each  segment  expand  each  factor  of  the  integrand  about  its  largest  term;  the  seg¬ 
ments  are  chosen  such  that  these  expansions  converge  as  rapidly  as  possible.  However,  this  is 
a  tedious  procedure,  and  one  which  must  be  adapted  to  the  relative  values  of  the  three  param¬ 
eters  involved,  namely,  r  /x  ,  L/\  ,  and  X  /x  .  Therefore,  in  the  present  work  we  shall  be 
,  j,  o'  e*  e  e'  m  ^ 

satisifed  with  the  approximation  (III-63). 

Integrals  such  as  (III-63),  of  the  form 


I(p,q) 


r 

n 


(1  +  V  )' 


dv 


(~1  <  p  <  2q  -  1) 


(III-64) 


are  easily  evaluated  in  terms  of  gamma  functions.  Substituting  v  =  tan</?,  (III-64)  becomes 
I(p,  q)  =  ^  (sin<^  (cos  d</? 


1-D/P+l  P+l\ 

=  2  'I- 

37 


where  B(x,  y)  is  the  beta  function  defined  by 

Therefore, 

1  r(P4i)r(q-^ 


Kp,  q)  = 

The  integral  in  (III-63)  is  then 


r(q) 


(111-65) 


and  inserting  this  result  in  (III-60)  we  find 
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(III-66) 


r^(T)  /Tr  \2 


1 6  nTtt 


\  A  / 


[The  gamma  function,  always  written  with  an  argument,  is  not  to  be  confused  with  F,  defined  by 
(111-27).] 

Using  (III-36)  and  (III-37),  and  evaluating  the  numerical  constant,  we  have  the  final  result 
for  the  dispersion 


6  =  0.145K  r  M 

loo  o 


(111-67) 


where  by  we  mean  the  rms  value 
For  a  TPF 

6  =  0.076A"^/® 


(III-68) 


on  an  rms  ripple  of  about  4®  in  zero  field.  This  is  certainly  a  reasonable  order  of  magnitude, 
although  it  must  be  remembered  that  the  value  chosen  for  (5  x  10  erg/cm  )  is  a  very  rough 
estimate.  It  should  be  kept  in  mind  that  this  dispersion  includes  not  merely  LMH,  but  compo¬ 
nents  with  wave  vectors  in  all  directions  and  of  all  wavelengths.  Therefore,  we  expect  it  to  be 
somewhat  larger  than  LMR  as  measured  by  Lorentz  microscopy. 

For  an  isotropic  film  (K  0),  K  A  i  M  H  in  (III-67);  as  the  external  field  H  0  the 
^  0  0*^0  ' 

dispersion  grows  very  large  until  (III-58)  is  violated.  This  is  a  long  wavelength  divergence 
(not  unlike  the  divergence  of  spontaneous  magnetization  in  a  two-dimensional  lattice  ).  The 
problem  here  is  that  the  Fourier  components  become  strongly  coupled  (see  Appendix  B),  and  the 
iteration  procedure  which  led  to  (III-39)  breaks  down.  In  fact,  for  an  isotropic  film  in  zero 
field  our  Fourier  transform  method  is  not  a  good  one,  and  an  approach  such  as  that  used  by 

Rother^^  is  better.  (His  main  result  is  a  ripple  with  mean  wavelength  proportional  to  r  and 

2.  ^ 

magnitude  proportional  to  K.r  /A,  where  the  constants  of  proportionality  are  ~10.)  Using 

^  ^  39  40 

another  method,  which  differs  considerably  from  ours,  Hoffmann  "  has  recently  obtained  the 
ripple  wavelength  given  by  (III-56)  and  amplitude  given  by  (III-67)  with  a  slightly  different  numeri¬ 
cal  factor. 

Summarizing  the  results  of  See.  III-E,  we  find  the  magnetization  dispersion:  (a)  to  be  typi¬ 
cally  a  few  degrees;  (b)  to  be  proportional  to  and  the  crystallite  size;  (e)  to  vary  as  the  — | 
power  of  the  uniform  anisotropy  and  of  h  ±  1  (for  fields  along  m^  in  the  easy  and  hard  directions, 
respectively);  and  (d)  to  depend  only  weakly  on  the  film  thickness. 
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IV.  UNIFORM  MAGNETIZATION  REVERSAL 


A.  INTRODUCTION 

Before  we  consider  the  dynamic  behavior  of  the  magnetization  ripple,  it  is  pertinent  to 
examine  the  process  of  uniform  magnetization  reversal,  i.e.„  coherent  rotation  with  M(r,t)  = 
m^(t) .  This  process  is  the  starting  point  in  the  treatment  of  nonuniform  reversal  to  follow,  and 
it  is  important  to  understand  how  it  occurs.  Furthermore,  earlier  treatments^*^’"^ ^  of  rotational 
switching  have  been  based  on  the  uniform  rotation  model  (but  with  a  phenomenological  damping 
term  included  in  place  of  detailed  knowledge  of  the  dissipative  processes  involved).  After  a  short 
description  of  the  uniform  rotational  mode,  we  will  formulate  the  equations  of  motion  in  an  aniso¬ 
tropic  film;  in  Sec.  V  they  will  be  used  in  obtaining  the  dynamic  response  of  the  ripple  (no  ex¬ 
plicit  solution  will  be  necessary  for  this  purpose).  Then  we  will  solve  these  equations,  with  and 
without  a  damping  term,  for  an  isotropic  film;  the  overdamped  case  will  be  treated  somewhat 
more  rigorously  than  has  been  done  previously. 

Rotational  switching  in  a  film  is  a  two-step  process.  When  a  pulse  magnetic  field  (typi¬ 
cally  a  few  oersteds)  is  applied  in  the  film  plane  at  time  t  =  0,  the  magnetization  m^(t),  initially 

lying  in  the  plane,  starts  to  process  about  H  .  In  doing  so  it  lifts  slightly  out  of  the  plane  (typi- 
-Z  ^ 

cally  10  radian),  creating  a  normal  demagnetizing  field  equal  and  opposite  to  the  normal  com- 
—  2 

ponent  of  "^^rn^  (typically  10  oe).  The  magnetization  now  processes  about  this  larger  field  in  a 
nearly  planar  path. 

However,  this  uniform  mode  has  acquired  an  energy  — H  •  from  the  external  field;  as 

a  result,  in  the  absence  of  any  dissipative  mechanisms,  m^  must  continue  to  process  indefinitely. 

If  there  is  a  slow  energy  loss,  m^  will  experience  damped  oscillations  about  its  equilibrium  direc- 


reaches  equilibrium  without  oscillations.  In  many  switching  experiments  oscillations  are  not  ob- 

42  10 

served  (but  see,  for  example,  Dietrich  and  Proebster  and  Hearn  ),  so  that  in  the  phenomeno¬ 
logical  theory  of  Sec.  IV-D  the  overdamped  case  is  of  great  practical  interest. 


B.  EQUATIONS  OF  MOTION 

The  gyromagnetic  equation  for  a  uniform  magnetic  moment  density  ni^(t)  is  [see  (1-2)] 
dm 

=  — ym^  X  h^  +  damping  term  (IV-1) 

where  y  is  the  absolute  value  of  the  gyromagnetic  ratio  and  h^  is  the  spatial  average  of  the 
effective  magnetic  field  in  the  ferromagnet.  For  a  thin  film  this  effective  field  may  be  taken 
as  the  external  field  (pulse  and  steady),  the  uniaxial  anisotropy  field,  and  the  uniform  demag¬ 
netizing  field —47rm^  •  (II-24a)]. 

As  long  as  we  remain  ignorant  of  the  specific  loss  mechanisms,  we  can  only  guess  at  the 

12 

form  of  the  damping  term.  Landau  and  Lifshitz  proposed 

- - — ^  m  X  (m  X  h  )  (IV-2) 

I  1 2  o  o  o 
m 

'  o ' 

where  \  is  a  positive  constant  of  dimensions  sec  A  slight  modification  was  suggested  by 
1 8 

Gilbert,  who  used  a  Lagrangian  approach  and  a  Rayleigh  dissipation  function  to  arrive  at 
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dm 
_ c 

dt 


dm 


=  — ym  X  h  + 
^  o  o 


m  X 
o 


dt 


(rv-3) 


where  a  is  a  dimensionless  positive  constant.  However,  it  is  easy  to  show  that  these  damping 
terms  are  equivalent,  except  for  a  contraction  of  the  time  scale.  Taking  the  dot  product  of  m^ 
with  (IV-3),  we  see  that 


dm 
_ c 

dt 


=  0 


(IV-4a) 


or 

I  m^l  =  const.  E  (IV-4b) 

Note  that  |rn^|  is  also  conserved  if  a  damping  term  (IV-2)  is  assumed  instead.  Now,  taking  the 
cross  product  of  m^  with  (IV-3),  using  (IV-4a,b),  and  rearranging  terms,  we  find 

y  dm 

(1  +  a  )  — =  -ym  Xh--^m  x(m  xh)  .  (IV-5) 

dt  '  o  o  M  o  o  o 
o 

2-1 

Thus,  if  we  let  t'  =  (1  +  a  )  t  in  (IV-5)  and  make  the  substitution 

\  =  ayM^  (IV-6) 


the  Landau-Lifshitz  and  Gilbert  forms  are  identical.  The  contraction  of  the  time  scale  is  unim¬ 
portant  if  a  «  1,  but  as  a  —  (IV-5)  gives 


whereas 


0 


(IV-7a) 


dm 

(IV-7h) 

In  other  words,  in  the  limit  of  large  damping,  the  switching  speed  goes  to  zero  for  Gilbert 

damping  but  becomes  infinite  for  LandaurLifshitz  damping.  The  latter  situation  seems  physically 
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unreasonable,  as  pointed  out  by  Kikuchi.  Callen  has  suggested  a  more  general  form  of  the 
damped  gyromagnetic  equation,  the  various  parameters  of  which  may  be  associated  with  specific 
physical  processes.  However,  the  difficulties  one  would  encounter  in  attempting  to  use  this 
approach  to  describe  a  s, witching  process  (in  which  m^  is  initially  far  from  equilibrium)  are 
probably  insurmountable.  The  Gilbert  equation  will  suffice  for  the  purpose  of  this  section,  and 
we  now  proceed  from  (IV-5), 

We  shall  use  a  spherical  coordinate  system  for  m^  as  shown  in  Fig.  7  where  the  +  x-direction 
is  the  easy  direction  of  uniaxial  anisotropy  nearest  m^(0) .  In  an  isotropic  film  we  take  rn^(O)  = 
M^i^.  If  we  let  of  rn^  in  the  external,  demagnetizing,  and  anisotropy 

fields,  then  the  effective  field  is 


m 

o 


E 

o 


(IV-8) 
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Fig.  7.  Caardinate  system  for  dynamic  case. 


and  the  spatial  components  of  (IV-5)  are 


^  de  aE  aE 

(1  +  — Q  3  _  y  o  _  ^  — o 

'  '  dt  Ttr  3<p  M  30 

o  ^o  o  o 


(IV-9a) 


and 


( 1  -h  a  )  sin  0 


d<^ 


aE 


aE 


O  _  y  _ o  _  Oiy _ o 


o  dt  M  ao  M  d(p 

O  O  0^0 


(IV -9b) 


We  consider  an  external  field  H  in  the  film  plane  at  an  angle  ^  to  the  x-axis.  Then  the  energy 
E^  is  given  by 

E  =  -m  •  H  +  27r(m  •  7  )^  -  K  M  "^(m  •  T 

o  o  oz  oo  ox 


2  2  2  2 

—  M  HsinG  cos(<^  —  Q)  ZnM  cos  0  —  K  sin  0  cos  </? 

o  o  o  oo  o  o 


(IV-10) 


with  derivatives 


1  1  2 
=  — H  cos  0  cos  iw  —  Q)  —  ZnM  sin  20  — ^  sin  20  cos  0 

M  a0  o  ^o  o  o  2  K  o 

o  o 


(IV-lla) 


1  1  2 

7^ — -  =  n  sin0  sin(g?  —  B)  +  ^  sin  0  sin  2(p 
M  a«>  o  ^o  ^  2  K  o  ^o 

o  ^  o 

where  the  anisotropy  field  is  defined  by 
2K 


H 


K  M 


(IV- 11b) 


(IV-12) 


We  next  let 


b  =  0  -  -y 

'^o  o  2 


(IV- 13) 


and  make  the  key  assumption  that  during  the  entire  switching  process  m^  is  confined  nearly  to 
the  film  plane,  or 
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(IV- 14) 


^  «  1 


It  is  easy  to  see  from  the  energy  expression  (IV-10)  that  the  condition  (IV- 14)  will  be  well  satisfied 
provided 

47rM  »H  (IV- 15) 

o 

because  the  magnetostatic  energy  created  in  lifting  m^  out  of  the  film  plane  an  appreciable  angle 
is  then  much  greater  than  the  energy  available  from  the  external  field.  The  requirement  (IV-15) 
is  indeed  met  for  all  experimental  situations  of  current  interest.  Also,  we  shall  verify  (IV- 14) 
for  the  solutions  to  be  found  in  Secs.  IV-C  and  IV-D.  The  condition  (IV- 14)  will  enable  us  to 
separate  the  coordinates  in  (IV-9a,  b),  but  even  if  it  were  violated  we  could  still  proceed  in  the 
next  section  with  the  calculation  of  the  dynamic  response  of  magnetization  ripple.  However, 
the  simplification  introduced  by  (IV- 14)  is  so  great  that  at  the  risk  of  some  loss  in  generality 
we  will  assume  it  holds  in  all  that  follows.  Then  with  the  further  requirement 

47rM  »  (IV-16) 

o  K 


which  is  invariably  met  in  practice,  (IV-lla,b)  become 


1 

M 

o 


dE 
_ o 

ae 

o 


47rM  lb 
0^0 


(IV- 17a) 


1 

M 

o 


^  H  sin{(p ^ 


sin  2,(0 
^  o 


We  next  generalize  (IV-17b).  Consider  a  step  function  pulse  field 


(IV- 17b) 


H  (t)  =  H  =  const.  t  >  0 
P  P 

=  0  t  <  0  (IV-18) 

at  an  angle  to  the  x-axis,  and  a  constant  bias  field  at  an  angle  Then  (IV- 17b)  becomes 

1  1 

for  t  >  0,  and  (IV- 17a)  is  of  course  unchanged.  We  note  at  this  point  that  in  the  usual  experimen¬ 
tal  arrangements  either 


or 

0  <  <  TT  ,  ^  ^  •  (IV-l9b) 

In  the  first  case  [(IV-19a)],  that  of  a  hard  axis  bias  field  and  an  easy  axis  switching  field,  the 
notation 


H  -  H|i 

P  II 

(IV-19C) 
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is  conventional.  In  the  second  case  (IV- 19b),  90®  switching  with  =  tt/Z  is  most  common.  Other 
arrangements  are  possible,  but  experimental  difficulties  invariably  render  them  undesirable. 
Inserting  (IV-17a)  and  (IV-17b-l)  in  (IV-9a,b),  and  using  (IV-13)  and  (IV-14),  wc  find 


(1  +  n  ) 


dt 


(IV-20a) 


where 


(1  +  a  ) 


d(p 
^  c 

dt 


u)  lb  +  aco.piw  ) 
m^o  h^  ^o 


(IV-20b) 


=  yHp  (IV-2  1a) 

a;  =  y47rM  (IV-2lb) 

m  o 

1 

-p(<Pq)  =  sin  ((P^  -  ^  )  +  ^  sin  ((p  -  ^  )  +  -2  sin  2(p^  (IV-21c) 

P  p  P  p 

The  precessional  frequency  in  the  pulse  field  is  and  is  typically  5  x  10^  sec”  cj^  is  the 

precessional  frequency  in  the  demagnetizing  field  for  m  perpendicular  to  the  film  plane  and  is 

11-1  ^  — 

1.9  X  10  sec  for  a  TPF;  pi(p is  proportional  to  the  planar  torque  on  m^.  It  is  convenient 

to  introduce  a  dimensionless  time  variable 


T  =  /o),  w  t  (IV-22) 

V  h  m 

and  to  normalize  the  angular  deviation  of  m^  from  the  film  plane  by  the  substitution 

a  =  e"%  (IV-23) 

0  ^0 

where 


£  =  /  —  «  1 
w  m 

The  small  parameter  c  is  typically  1.5  X  10 
parameter 


a 

£ 


Equations  (IV-20a,b)  then  become 


-2 


2  2 

(1  +  €  V  )  =  p((p^)  -  va^ 


2  2.  2 
(1  +  e  V  )  (pQ  =  +  e  vp(ip ^ 


(IV-24) 

We  also  define  a  normalized  damping 

(IV-25) 

(IV-26a) 

(IV-26b) 


where  the  dot  indicates  differentiation  with  respect  to  t.  For  the  purposes  of  Sec.  V,  the  phe¬ 
nomenological  damping  will  be  small  {p  <  1),  and  (IV-26a,b)  may  be  approximated  by 


(y  =  p{(p  )  —  VO 
o  ^  ^  o  o 


=  a 

o 


(IV-27a) 

(IV-27b) 
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The  qualitative  description  of  rotational  switching  given  in  Sec.  IV-A  is  reflected  in  {IV-26a,b) 
[or  (IV-27a,  b)].  The  primary  precession  about  the  external  field  (plus  uniaxial  anisotropy  field) 
is  described  by  (IV-26a),  and  the  secondary  precession  about  the  normal  demagnetizing  field  by 
(IV-26b).  We  see  that  it  is  mainly  the  primary  precession  which  is  damped;  this  leads  to  a  re¬ 
duced  demagnetizing  field  and  therefore  a  slower  secondary  precession. 

Initial  values  of  and  (p ^  are  determined  by  minimizing  the  energy  for  t  <  0.  Equations 
(IV-17a)  and  (IV-17b-l)  give  [using  (IV-23)] 

CJ^(O)  =  0  (IV-28a) 


sin  [<^^(0)  -  /3^]  +  2  sin  2(p ^{0)  =  0 
with  the  added  condition 


i  O 


A  2 

O  d(p 
^  o 


cos  [<^^(0)  -  /3^]  +  cos  2.(p^{0)  >  0 


t<0 


For  the  usual  arrangement  =  7r/2,  (IV-28b)  and  (IV-29)  give 

H 

(P^{0)  =  arc  sin  yj-  (0  ^ 


K 


K' 


(Hb^HK) 


(IV-28b) 


{IV-29) 


(IV-30a) 

{IV-30b) 


[The  case  (IV-30b)  is  trivial,  as  no  switching  occurs.] 

To  conclude  Sec.  IV-B,  we  examine  the  behavior  of  the  uniform  mode  inthe  vicinity  of  its  new 
equilibrium  direction  at  the  completion  of  the  switching  process.  This  direction  is  determined  by 


lim 

t-^OO 


=  0 


{IV-31) 


or,  equivalently,  that  the  energy  is  minimized  for  and  (p ^  at  t 

lim  CT  (t)  =  cj  =0 
.  o  o 


oo.  Either  way,  we  find 
{IV-32a) 


lim  p[(^  (t)]  E  p(^  )  =  0  .  {IV-32b) 

t-^OO 


We  expand  (IV-26a,b)  about  expand  and,  assuming  that  (J  ^  and  ~  <^q)  vary  as 
exp[fi(l  +  ^t],  arrive  at  the  characteristic  determinant 


~P  * 

—  1  —  6  ^i/p 


=  0 


{IV- 3  3) 


where 


^0=^0 


Dropping  terms  of  order  e  ,  we  extract  from  {IV-33)  the  critical  damping 


{IV-34) 
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(IV-35) 


V  =2  V  — p' 
c  ^ 


or,  with  the  help  of  (IV-21c)  and  (IV-24) 

2 


UttM 
V  o 


[H  cos  {(P^  -  p  cos  cos  2(p  ^ 


,1/2 


(IV-36) 


If  a  <  a  ,  m  overshoots  equilibrium  and  experiences  a  damped  oscillation;  if  cv  >  o  ,  m 

^  V 

reaches  equilibrium  without  oscillations.  Note  that  if  p'  >  0,  the  equilibrium  direction  0^  is 
unstable;  thus  the  threshold  field  for  irreversible  rotation  through  any  equilibrium  point  is 
determined  by 


p  =  p 


=  0 


(IV-37) 


(See  p.  22.) 

C.  UNDAMPED  UNIFORM  MODE 

An  independent  measurement  of  a  may  be  obtained  from  the  linewidth  in  ferromagnetic 

resonance.  Although  there  may  be  dissipative  processes  in  large-angle  switching  which  are  not 

operative  in  small-angle  resonance  (and  indeed  it  is  one  purpose  of  this  work  to  uncover  such 

processes),  the  microwave  resonance  damping  constant  ^ provides  a  measure  of  "intrinsic" 

loss  mechanisms  (the  details  of  which  we  will  make  no  attempt  to  treat)  and  a  lower  bound  on  the 
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total  damping.  Various  measurements  in  Permalloy  films  at  microwave  frequencies  '  ' 

_2 

give  a  isi  10  ,  which  is  barely  less  than  a  for  the  usual  switching  experiments. 

As  an  approximation  to  this  slightly  underdamped  situation,  we  solve  (IV-27a,  b)  with  p  =  0 
for  an  isotropic  film  without  a  bias  field.  The  equations  of  motion,  from  (IV-21c)  and  (IV-27a,  b), 
are 


a  =  —sin  {w  —  (3) 
o 


0  =  G 

^  o  o 


(IV-38a) 

(IV-38b) 


with  initial  conditions 


(^^{0)  =  0^{O)  =  0 


(rv-38c) 


Eliminating  from  (IV-38a,  b),  we  obtain 

0  +  sin  {0  —  ^)  =  0 

^  o  ^  o 


(IV-39) 


with  initial  conditions 


<P^iO)  =  0^{O)  =  0 


(IV-40) 


Equation  (IV-39)  is  the  equation  of  a  simple  pendulum,  and  is  readily  solved  by  Jacobian  elliptic 
44 

functions.  We  introduce  a  new  dependent  variable 


u  =  b  tan  2  iP  -  (P^) 


(IV-41) 


where 


b  =  cot  - 


{IV-42) 
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Equations  {IV-39)  and  (IV-40)  then  become 


2uu^  1  A 
u--^^  +u  =  0 

b  +  u 


(IV-43a) 


u{0)  -  1  (rV-43b) 

u(0)  -  0  .  (IV-43c) 

An  integrating  factor  for  {IV-43a)  is  (b^  +  u^)  from  which  we  find  the  first  integral 


^2  ^  (b^  4-  u^)  (1 
b^  4-  1 


(IV-44) 


A  second  integration  gives  us 


u  =  cn{T,  sin-^) 


{rV-45) 


which  is  the  main  result  of  Sec.  IV-C. 

As  expected  for  a  conservative  system,  the  motion  is  periodic,  with  period  4K(sin/?/2), 
where  K  is  the  complete  elliptic  integral  of  the  first  kind.  The  switching  time  may  be  de¬ 
fined  as  the  time  it  takes  m^  to  rotate  from  =  0  to  (p ^  =  (3,  the  new  equilibrium  angle,  and  is 


T  =  K(sin 
s  2 


{rv-46) 


K(sin|-) 


(IV-47) 


Finally,  we  must  verify  (IV-14).  From  (IV-38b)  and  (IV-41), 
-2b 


a  =  ^ 

o  ,  2 


b  +  u 

and  using  (IV -4 5),  we  find 


([V-48) 


O'  =  sinfi 
o  ^ 


i  • 

sn{T,  sin-^) 
dn(T,  sin  y) 


{IV-49) 


This  function  attains  its  maximum  value  at  r  =  [see  {IV-46)],  and  is 


a  (t  )  =  2  sin  4 
o  s  2 


{IV-50) 


Thus 


{rv-51) 


D.  OVERDAMPED  UNIFORM  MODE 

In  contrast  to  the  damping  deduced  from  microwave  resonance  linewidths,  the  damping 
actually  observed  in  rotational  switching  experiments  is  quite  high,  particularly  in  the 
intermediate-speed  region  (see  p.3),  where  one  invariably  finds 


38 


V  »  1 


(IV-52) 


f«2» 


H 

_ e_ 

47rM 

o 


We  leave  for  Sec.  VI  the  question  of  why  this  damping  should  be  so  great  and  solve  (IV-26  a,b), 
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for  an  isotropic  film  with  a  damping  constant  satisfying  (IV-52),  by  a  boundary  layer  method. 
We  first  eliminate  from  (IV-26a,b),  obtaining 

{i  +  +  v<p^-  p  =  0  (IV-53) 


where  for  an  isotropic  film 


p(<Po)  =  -sin(<p^  -  /3) 


^q(O)  =  0 


*  /n\  _  f  ^  sin^ 

1  +  e  r 


Next  we  introduce  a  new  time  variable 


(IV-54a) 
(IV- 54b) 

(IV-54C) 


-1 

T  =  V  T 


i^it 


a 


(IV-55) 


and  transform  (IV-53)  to  an  equation  for  u(t),  where  u  is  defined  by  (IV-41).  The  result  is 


6 


2u  .  du  .  2 

,2  ^  2  ^df  ^ 

b  +  u 


du 

dr 


+  +  u  =  0 


(IV-56a) 


with 


and 


where 


u(0)  =  1 


du 

dr 


- 

=  —  e  6 


T=0 


1  +  a 


(IV- 56b) 


(IV-56c) 


6  = 


^  X  2  2 

1  +  e  r 


=  -^  (1  +  o^)  «  1 
a 


(IV-57) 


We  start  our  analysis  with  (IV-56a-c).  It  might  seem  at  first  glance  that  since  6  is  a 
small  parameter,  we  could  neglect  the  first  term  of  (IV-56a),  retaining  the  second  (damping) 
and  third  (restoring  torque)  terms.  This  is  the  "viscous  flow"  approximation  of  Smith"^  (which 
reduces  the  problem  to  an  integration  for  the  anisotropic  case).  However,  this  solution  is  not 
valid  in  the  neighborhood  of  t  =  0  because  there  (unless  a  >  1)  du/df  is  small  and  the  first 
term  in  (IV-56a)  is  of  order  unity.  This  boundary  layer  at  t  =  0  must  be  accounted  for  properly 
in  order  to  apply  the  initial  conditions  (IV-56b)  and  (IV-56c)  to  the  "viscous  flow"  solution  valid 
past  the  boundary  layer. 

Let  the  solution  beyond  the  (as  yet  unspecified)  boundary  layer  be  u(t),  which  we  expand  as 
a  power  series  in  6: 
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oo 


u(T)  =  Yj  ■5'  (IV-58) 

i=0 

Inserting  (IV-58)  in  (IV-56a),  we  find  that  the  solution  to  zeroth  order  in  6  is 

u  =  C  (IV-59) 

o  o 

where  is  a  constant  which  must  be  determined.  Higher  order  terms  u^,u^,  .  .  .  can  be  found 
without  difficulty,  but  we  shall  not  give  them  explicitly.  It  is  obvious  from  {IV-56a)  that  they  will 
each  contain  a  single  undetermined  constant  C^,  .  .  .  . 

To  solve  (IV-56a-c)  within  the  boundary  layer,  we  first  expand  the  time  scale  through  the 
change  of  variable 

V  =  (m  <  0)  {IV-60) 

where  m  will  be  chosen  such  that  the  boundary  layer,  in  which  the  first  term  of  {IV-56a)  is  not 
negligible,  is  given  by  0  ^  v  S  1.  We  next  let  the  complete  solution  to  {IV-56a-c),  valid  for  all 
time,  be 


u  =  u(t)  +  6  w(v)  (n  >  0)  {IV-61) 

where  n  will  be  chosen  such  that  w  =  0(1)  for  0  <  v  ^  1.  Since  u  -*  u  for  v  »  1,  we  require 


lim  w(v)  =  0 

V  ^  oo 


(IV-62) 


With  these  new  variables,  (IV-56a-c)  become 

,  ,.2m+n+l  ,,  (u  +  6*^w)  (u'  +  6^^^w’)^  ,  ,  .m+n  ,  ,  ^  ,  ^n  ^  , 

6u"  +  6  w* 26  ^ — — - R - —  +  u'  +  6  w'  +  u  +  6w-0  (IV -6  3a) 

b  +  (u  +  6*^w) 


u(0)  =  u(0)  +  6'^w(0)  =  1  (IV-63b) 

u-(0)  =  u-(0)  +  >  -  1  (IV-63c) 

1  +  a 

where  u'  =  du/dr,  and  w’  =  dw/dv.  Now  let  us  assume  that  w  is  a  well-behaved  function  so  that 
w'  and  w”  are  of  order  unity  for  0  <  v  <  1.  Then  the  boundary-layer  width  is  determined  by  the 
condition  that  coefficients  of  w'  and  w”  in  (IV-63a)  are  of  equal  order  in  6,  giving 

m=-l  (IV-64) 

The  initial  values  of  w  and  w'  are  found  from  (IV-58),  (IV-61),  (IV-63b),  and  (IV-63c)  to  be 


w(0)  =  6"'^  [1  -  u^(0)  -  6u^(0)  -  .  .  .  ]  (IV-65a) 

and 

W'(0)  =  -  U^(0)  -  6u^(0)  -  .  .  .  ]  .  (IV-65b) 

The  exponent  n  and  the  constant  of  (IV-59)  are  then  determined  through  the  requirement  that 
w(0)  and  w’(0)  remain  of  order  unity  as  6  —  0.  From  the  first  two  terms  in  (IV-65a)  we  must 

have  either  n  =  0  or  u  (0)  =  1.  However,  if  n  =  0,  the  leading  terms  in  w'(0)  are  of  order  6  and 

2  ^  — 

€  ,  in  violation  of  the  above  requirement.  Therefore  u  (0)  =  1,  or  from  (IV-59), 
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1 


(IV-66) 


But  now  u^(0)  =  —  1,  so  that  the  largest  term  in  (IV-65b)  is  6^  ^  (remembering  that  6  <  1), 

from  which  we  see  that 

n  =  1  .  (IV-67) 

Next  we  expand  w  as  a  power  series  in  6: 


/(v)  =  Yj  w.(v)  6^ 
i=0 


(IV-68) 


Returning  to  (IV-63a)  with  m  =  —  1  and  n  =  1,  we  find  to  zeroth  order  in  6 


w* *  4-  w ’  =  — u*  —  u  =0 
o  o  o  o 


(IV-69a) 


with  initial  conditions  given  by  (IV-65a,b): 


w^(0)  =  -u^(0) 


w '  (0)  =  — u’  (0) 
o  o 


2  2 
1  +  q:  1  +  a 


(IV-69b) 

(IV-69C) 


The  solution  to  (IV-69a)  is 


w  =  D  e  +  D 
o  o  o 


(IV-70) 


where  and  are  constants  determined  by  (IV-69c)  and  (IV-62),  respectively.  The  result  is 


...  2-1  -V 

w  =  —(1  +  O'  )  e 
o 


(IV-71) 


The  condition  (IV-69b)  determines  u^(0),  i.e.,  the  constant  C^.  Then  from  the  next  term  in 
(IV-68),  comes  the  constant  in  u^,  and  so  on  to  any  desired  order. 


^2  2^ 

Collecting  our  results,  wc  have  found 
u  =  e"^  +  0(6) 


(IV-72a) 


and 


du  -T  , 

1 —  —  e  4"  ~ 

i  + 


- -  e  +  0(6) 


(IV-72b) 


where  the  second  term  in  (IV-72b)  is  of  importance  only  within  the  boundary  layer  defined  by 


0  T  £  6  <<  1 
or  in  real  time  by 


(IV-73a) 


0  ^  t  < 


1  -I-  g 
y47rM^g 


(IV-73b) 


The  switching  time  for  this  overdamped  case  may  be  defined  conveniently  as  the  decay  time  of 
u,  which  means 


T  =  1 

s 


(IV- 74a) 
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or  in  real  time 


t  -  — 
s  yH 

(IV-74b) 

From  (IV-26b),  (IV-41),  (IV-55), 

and  (IV-72a,  b) 

1  2b 

(e-"  -  . 

(IV-75) 

^o  ^  ,  -2t 

b  +  e 

Thus  from  (IV-23)  and  (IV-42), 

^^o^max  ~ 

for  n  >  p  >  ^ 

-  1  sinp 

for  ^  ^  /S  >  0 

(IV-76) 

and  (IV-14)  is  satisfied.  A  simple  physical  interpretation  can  be  found  for  the  boundary  layer: 
it  is  the  initial  interval  during  which  the  growth  of  the  component  of  m^  normal  to  the  film  plane 
takes  place,  i.e.,  the  time  for  the  primary  precession.  [From  (IV-73b)  we  see  that  this  time 
has  a  minimum  at  a  =  1  of  z/ w  ~  10  ^  ^  sec.l 

The  qualitative  results  we  have  obtained  here  for  overdamped  switching  (and  in  Sec.  IV-C 
for  undamped  switching)  will  still  be  valid  for  anisotropic  films.  But  if  accurate  numerical 
results  are  desired,  the  presence  of  anisotropy  (except  for  some  limiting  cases)  so  complicates 
the  integrations  that  a  computer  solution  becomes  advisable.  However,  a  rather  good  approx¬ 
imation  is  to  replace  H  in  the  solutions  for  </?  and  a  by  H  —  H  where  H  ^  is  the  threshold 
field  for  irreversible  rotation  [see  (IV-37)]. 
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V.  SPIN-WAVE  THEORY:  QUASISTATIC  SPECTRUM 
AND  TRANSIENT  RESPONSE 


A.  INTRODUCTION 

We  now  turn  to  the  dynamic  behavior  of  magnetization  fluctuations,  i.e.,  spin  waves,  paying 
particular  attention  to  the  transient  response  of  ripple  to  a  pulse  field;  we  also  obtain  the  quasi¬ 
static  spin-wave  spectrum  and  briefly  discuss  instabilities.  The  microwave  resonance  situa¬ 
tion,  in  which  the  amplitudes  of  uniform  and  nonuniform  modes  are  small,  is  rather  well  under¬ 
stood  in  the  short-wavelength  limit^^'^^  (in  which  magnetostatic  fields  for  all  modes  with  k  0 

may  be  replaced  by  their  infinite  medium  values)  and  in  the  long-wavelength,  or  magnetostatic 
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limit  '  (in  which  exchange  fields  may  be  neglected).  Our  treatment  differs  from  these  in  two 

important  respects.  First,  the  amplitude  of  the  uniform  mode  cannot  be  considered  small  (it 

may  be  180°)  and  must  therefore  be  retained  to  all  orders.  We  avoid  some  of  the  difficulty  this 

nonlinearity  introduces  by  using  a  coordinate  system  rotating  with  the  uniform  mode.  Second, 

neither  the  infinite  medium  nor  the  magnetostatic  approximations  are  adequate  for  a  thin  film 

in  low  fields,  since  the  spin  waves  most  important  in  switching  have  wavelengths  somewhat 

greater  than  the  film  thickness  and  also  involve  exchange  fields.  In  resonance  also,  these  ap- 
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proximations  break  down.  This  difficulty  is  surmounted  with  the  help  of  the  thin-film  approxi¬ 
mation  (TFA)  of  Sec.  II,  which  enables  us  to  consider  magnetostatic  and  exchange  fields  simul¬ 
taneously.  (This  has  sdready  been  used  to  obtain  the  static  solution  of  Sec.  III.) 

We  shall  assume  that  random,  local  anisotropy  forces  have  no  first-order  dynamic  effects. 
Since  these  forces  are,  on  the  average,  isotropic,  their  dynamic  perturbing  effects  should  be 
approximately  as  great  as  their  static  effects;  but  we  shall  find  that  magnetostatic  fields  grow 
very  large  during  magnetization  reversal,  so  that  we  may  expect  local  anisotropy  fields  to  be 
relatively  unimportant.  This,  together  with  the  assumption  that  fluctuations  of  M  from  m^  are 
small  [see  (I-la)],  implies  that  ripple  components  are  uncoupled  dynamically  (as  well  as  stati¬ 
cally),  since  the  only  linear  coupling  is  via  spatially  varying  anisotropy  fields,  as  was  shown  in 
Sec.  III. 

Another  assumption  we  make,  in  order  to  avoid  excessive  complications,  is  that  damping 
may  be  neglected.  As  far  as  the  uniform  mode  is  concerned,  the  intrinsic  (resonance)  damping 
is  small,  as  was  noted  in  Sec.  IV;  the  large  damping  observed  in  switching  is  precisely  what  we 
are  attempting  to  find  a  physical  basis  for,  and  we  therefore  avoid  any  phenomenological  treat¬ 
ment  of  it.  Damping  of  spin  waves  can  safely  be  neglected  if  it  results  in  relaxation  times  which 
are  longer  than  the  switching  time.  If,  on  the  other  hand,  relaxation  times  are  short,  the  initial 

Fourier  components  do  not  maintain  their  identity  during  switching,  and  the  problem  is  quite 

-4 

different.  Since  the  important  components  are  relatively  long  wavelength  10  cm),  we  may 
expect  their  damping  to  be  about  the  same  as  that  of  the  uniform  mode  in  resonance,  and  we  are, 
therefore,  justified  in  neglecting  it. 

Finally,  we  assume  that  the  z-dependence  of  6M(r,t)  may  be  ignored  (z-axis  perpendicular 

to  the  film  plane).  We  saw  in  Sec.  Ill  that  exchange  fields  sharply  attenuate  all  components  with 

wavelengths  shorter  than  an  exchcinge  wavelength  [see  (III-36)]  which  is  typically  an  order 

of  magnitude  greater  than  the  film  thickness.  Thus,  any  components  with  appreciable  z- 

dependence  would  be  expected  to  have  very  small  amplitudes.  In  somewhat  thicker  films 
o  50 

(^2000  A),  standing,  z-directed  spin  waves,  which  can  be  excited  by  microwave  fields  and  which 
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depend  on  the  boundary  conditions  for  6M  at  the  film  surfaces,  might  play  a  role  in  magnetization 
reversal;  we  will  not  consider  this  possibility  in  the  present  work.  (Note  also  that  in  our  model 
for  the  microstructure  of  a  polycrystalline  film,  local  anisotropy  forces  are  constant  across 
the  film  thickness,  precluding  z-dependent  dispersion-induced  fluctuations  of  M.) 

B.  EQUATIONS  OF  MOTION 

The  gyromagnetic  equation  for  the  magnetization  M(r,t)  of  a  ferromagnet,  as  given  by  (1-2) 
without  damping,  is 


aM(r,  t) 

at 


=  -YTgff(r,  t) 


(V-1) 


where 


f  -  M(r,  t)  X  t)  .  ( V-2) 

Separating  the  spatially  varying  part  of  M  by  (I- lb)  and  by 

ileff  =  h^ii)  +  <5H(?,t)  (V-3) 

(where  b^(t)  is  the  spatial  average  of  H^^^(r,t)J,  and  taking  a  spatial  average  of  (V-1),  we  find 
dm  _  _  _ 

-^  =  -Ym^  X  -  y<6M  X  6H>  (V-4a) 

This  equation  for  the  uniform  mode,  with  the  second  term  on  the  right  side  replaced  by  a  phe¬ 
nomenological  damping  term,  was  discussed  in  Sec.  IV.  The  nonlinear  reaction  of  spin  waves 
on  the  uniform  mode  is  described  by  this  second  term  and  will  be  computed  in  Sec.  VI.  Sub¬ 
tracting  (V-4a)  from  (V-1),  we  obtain  the  dynamic  equation  for  fluctuations  of  the  magnetization 

6M  =  -Y(m^  X  6H  -  X  6M)  -  y(6M  X  6H  -  <6M  x  6H»  (V-4b) 

Since  (as  we  shall  find)  6H  is  proportional  to  6M  (or,  strictly  speaking,  to  an  integral  operator 
on  6M),  the  assumption  (I- la)  enables  us  to  neglect,  to  a  first  approximation,  the  second  term 
on  the  right  side  of  (V-4b),  which  contains  the  effects  of  interactions  among  the  components  of 
6M. 

This  becomes  clearer  when  we  expand  6M  and  6H  in  Fourier  series: 


6M(r,  t)=  Yj 

k^O  ^ 

6H(r,t)  =  Y  • 

k^O  ^ 

Then  (V-4a)  becomes 

dn^o  -  -  V  -  - 

=  — ym  X  h  —  y  /,  m_^  ^  h  _ 

dt  ^  o  o  ^  r 


(V-5a) 


(V-5b) 


(V-6a) 
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and  the  k  ^  0  components  of  (V-4b)  are 
dm 


=  -y{m  Xh^-h  X  m_^)  -y  ^  m_^  ^  ^  h 

^  k-  ^  k-  k-k-'  k* 

k'7^k,0  ^  ^ 


(V-6b) 


We  will  show  below  that  the  k^^  component  of  the  field  is  related  to  the  k^^  component  of  the 

magnetization  by  a  field  dyadic  K_^ 

k 


m_^(t) 

h  (t)=F  (t)- 
k  k  o 


(V.7) 


Then  the  same  iterative  procedure  that  was  used  to  solve  (III-Z3a,  b)  may  be  applied  to  (V-6a,b), 
with  the  resultant  first  approximations 


dm 

=  — ym  X  h 
dt  o  o 


(V-8a) 


and 


dm 


dt 


=  — y(m  X  h  —  h  x  rn  ) 
°  k  °  k 


(V-8b) 


We  use  a  coordinate  system  based  on  m^,  as  shown  in  Fig.  7,  in  which  the  spatial  compo¬ 
nents  of  m  are  written 
k 


m  =  M  ( i  m_  +  i^0  +  i  cp  )  =  M  (m  ,0  ,  (p  ) 

k  °  k  ®  k  k  °  k  k  k 


(V-9) 


With  the  spatial  components  of  h_^  given  by 

k 


(all  k) 


k  k  k  k 
and  using  the  relations 


(V-10) 


i  =  {0,0  ,<p  sin  0  ) 
m  '  o  o  o 


i^  =  (—0  ,  0,  <p  cos  0  ) 
0  '  o  o  o 


(V-lla) 

(V-llb) 


=  {—^  sin0  ,—(p  cos  0  ,  0) 
(p  ^  o  o' 


(V-llc) 


(which  are  most  easily  found  by  inspection  of  Fig.  7),  we  obtain  for  the  i^-component  of  (V-8b) 


dm 


dO  d<p  ^ 

-TT —  =  — 0  +  “-rp  sinO  <p  —  y(h^O  —  h  ^(p  ) 

dt  dt  dt  o 


(V-12) 


t  If  spatially  varying  farces  such  as  local  anisatrapy  are  present,  (V-7)  has  an  additional  term  M  I  ' 

and  the  camponents  are  coupled.  k*  k-k’  k' 
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However,  (V-8a),  in  terms  of  the  field  components  (V-10),  becomes 


d0 


dt 


°  =  Yh*^ 
o 


(V-13a) 


d^ 

dt 


o  .  ^  ,0 

sin  O  =  — vh 
o  o 


and  inserting  these  in  (V-12)  we  find 
dm 

=  0 


dt 


(V-13b) 


(V-14) 


Since  m  is  initially  zero, 
k 

m  (t)  -  0 
k 


(V-15) 


in  our  first  (linear)  approximation.  [We  could  also  have  obtained  (V-15)  from  the  constraint 
|m|  =  const.  =  M^,  which  follows  from  (V-1)  and  (V-2).)  From  (V-13a,  b)  and  (V-15),  the  and 
i^  components  of  (V-8b)  arc  then 


d0 

_ [ 

dt 

dep 


^  =  — y(h  ^  +  h  ^  cot  0  )  <^  +  yh  ^ 

^  k  k 


dt 


-  y(h  +  h  cot  0)0  —  yh 

«  k  k 


(V-I6a) 


(V-l6b) 


The  exchange  contribution  to  h  is  found  from  (1-4) 

k 


H  =  V^M 
o 


with  Fourier  components 

'  k'e  M 


(V-17) 


(V-18) 


The  magnetostatic  contribution  is  given  in  the  TFA  by  (II-24b): 


(hj  =-4ir\ 
'  k'm  I 


kk 


=  -4jr|^  X  (kL)  +  i^i^x(kL)|  •  m 

k 


where 


x(k)  =  K  e  ^  sinh  K  (=  1  —  k  4.  .  .  for  k  «  1) 
x(k)  =  1  -  x('^) 

Finally,  the  effective  uniaxial  anisotropy  field  [see  (III-2)]  is 


(=  K  —  K ^  .  for  K  «  1 ) 


(V-19a) 

(V-19b) 

(V-19c) 


M 


H  =  IL,  i  i  T,- 

a  K  X  X  M 


(V-20) 
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with  Fourier  components 


(h  )  =  i  i 

'  o'a  K  X  X 


(V-21a) 


H^i  i 
K  X  X 


(V-21b) 


We  see  from  (V-18),  (V-19a),  and  (V-21b)  that  h_^  can  indeed  be  expressed  in  the  form  (V-7), 
with  ^ 


K 

k 


-47rM 


2A 

M 


I  +  i  i 

K  X  X 


(V-22) 


where  I  is  the  idemfactor,  or  unit  dyadic.  The  vectors  in  (V-22),  in  the  (m,0,</?)  coordinate 
system,  are 


k  — ^ 

^  =  4(4»)  =  [sinG^  cos[(p^  -  4»),  cosG^  cos  -  4*),  -sin  -  4*)]  (V-23a) 

i^  =  (cosG^,  — sinG^,  0)  (V-23b) 

=  T(0)  (V-23c) 

where  the  wave  vector  k  lies  in  the  film  plane  at  an  angle  4*  to  the  x-axis.  We  also  note  that 
(V-l6a,b)  contain  the  factor 

h  ^  +  h  ^  cot  G  =cscG  h^  (V-24) 

o  o  o  o  o 

where  h^  is  the  component  of  the  uniform  field  along  the  projection  of  m^  onto  the  film  plane. 

It  is  convenient  to  write  (V-l6a,b)  in  the  form 

dG^ 
k 

-77 -  =  —00  ^G  —  CjO 

dt  (pO  -  k 

d(p 

V  V 

k  k 


( V-25a) 


( V-25b) 


For  the  case  of  an  external  field  H  in  the  film  plane  at  an  angle  /3  to  the  x-axis,  and  m^  nearly 
in  the  plane  such  that 

we  find,  with  the  help  of  (V-2  la) ,  (V-22),  (V-23a-c),  and  (V-24), 

2  2  2  ~  2  2 

<^00  =  cos(<p^ -/3)  +  cos  +  WgL  k  +  cos  ((^>q-4>)+x]  (V-27a) 

2  2  ~  2 

00  =00.  cos  {(p  —  /3)  +  oj  cos  2(p  +a)Lk  +co  x  sin  (ip  —  4>)  (V-2  7b) 

ipip  h  '^o  a  ^o  e  m'^  '^o  '  '  ' 
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(V-27c) 


=  CJ  ^ 
0(p  (pO 


1 

•=•  ip  [—a)  sin2(/?  +  a)  y  sin2((/?  —  4>)] 


where 


=  y|h| 

( V-28a) 

(V-28b) 

2A 

"^e  "  2 

(V-28c) 

ML 

o 

io  =  y47rM 
m  o 

(V-28d) 

Note  that  our  definition  of  the  exchange  frequency  is  not  the  conventional  one,  since  it  is  based 
on  the  half-thickness  of  the  film  L,  rather  than  on  the  lattice  constant  f  .  Equations  (V-13a,  b), 
which  describe  the  uniform  mode,  now  become 


1 

=-‘"h  sin(<p^-^)-co^  2 
d(p 


dt 


^  ^  f^u  COS((p  —  i3)  +  CJ  COS^  (p  +  <jO  ] 

o  ^  h  '^o  a  ^o  m^ 


(V-29a) 


(V-29b) 


The  coupled,  linear  equations  (V-25a,b),  in  which  the  coefficients  are,  in  general,  time 

dependent  through  the  time  dependence  of  (p^  and  ip^,  describe  the  motion  of  the  components  of 

6M  out  of  the  film  plane  (©_)  and  in  the  plane  {(p  ).  In  Secs.  V-C  and  V-D  we  will  find  quasistatic 

k  k 

spin-wave  solutions  (with  jp^,  (p ^  constant)  and  transient  solutions  (with  H  turned  on  at  t  -  0). 

C.  QUASISTATIC  SPECTRUM 

Here,  we  assume  that  there  exists  a  time  interval  over  which  di  (t)  and  (p  (t)  may  be  con- 

o  o 

sidered  constants  while  0_^(t)  and  <^_^(t)  are  rapidly  varying.  Then  solutions  of  (V-25a,  b)  are 

k  k 

spin  waves  with  time  dependence 


ict)  t 
k 


(V-30) 


where  the  eigenfrequencies  gj  are  given  by 

k 


2  2 
QO  (p(p  0(p 


(V-31) 


We  first  examine  the  spectrum  (V-31)  for  the  usual  resonance  situation  of  oscillations  about 

V 

a  stable  equilibrium.  Then  ip^  and  (p^  are  given  by  their  equilibrium  values  ip^  and  (p^  [which 
may  be  found  from  (V-29a,b)] 

^  =  0 


H  sin  {(p^  -  /3)  +  -  sin  Zcp^  =  0 


(V-32a) 

(V-32b) 
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with  the  stability  condition 


H  cos  (<^^  -  ^)  +  Hj.  cos  2^^  >  0  (V-32c) 

V  V 

Inserting  and  (p^  in  (V-27a-c),  (V-31)  becomes 

2  V  2  V  2  2 

^  cos(<^  —  P)  +  ct)  cos  (D  fcoLk  +0J  y] 

^  h  '^o  a  ^o  e  m'^  ^ 

k 

X  [W|^  cos(^^ -/3)  +  cos2^^  +  +  W^X  sin^(^^  -  4>)]  (V-33a) 

The  eigenfrequency  of  the  uniform  mode  may  be  found  by  expanding  (V-29a,b)  about 

(p^  =  or  from  (V-33a,b)  in  the  limit  k  -**  0,  and  is  given  by 

(jj^  -  fo),  cos  ((S  —  p)  +  u)  cos^  w  +  <jo  ]  [w,  cos  (w  —  p)  +  (jj  cos  Zw  1  ( V-33b) 

o  ^  h  '^o  a  ^o  m^  ^  h  '^o  a  ^o^  '  ' 

If  the  external  field  is  much  greater  than  the  anisotropy  field,  =  p,  and  (V-33a,b)  become 
(for  H  along  i^) 

(V-34a) 

(V-34b) 


2  2  2  2  2  ~  2 
^  +  w  L  k  +  cj  y)  (cj,  +  cj  L  k  +  a;  y  sin 

'  h  e  m'^'  '  h  e  m^  ' 

k 

o  h'  h  m' 


which  are  shown  in  Fig.  8  for  a  TPP^  in  a  field  of  800  oe.  In  the  infinite  medium  limit  kL  —  *5 

OJ  ^  -  (w.  +  w  L^k^)  (w,  +  CJ  iJk^  +  w  sin^  4>)  ( V-35) 

r*“  he  he  m  '  '' 

k 
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which  is  the  dispersion  relation  for  a  thin  disk  magnetized  in  its  plane  as  given  by  Callen.  Ex¬ 
tending  this  spectrum  to  k  =  0  (see  Fig.  8),  we  see  that  the  uniform  mode  lies  at  the  top  of  the 

spin-wave  manifold  (4>  =  tt/Z),  in  contrast  to  the  TP"A  spectrum  (V-34a),  which  collapses  to  a 
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point  as  k  —  0.  We  also  show  in  Fig.  8  the  magnetostatic  modes  found  by  Damon  and  Eshbach. 

Since  the  location  of  the  uniform  mode  relative  to  the  spin-wave  manifold  determines  the  coupling 
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of  the  uniform  mode  to  spin  waves,  '  the  distortion  of  the  infinite  medium  spectrum  at  long 

wavelengths  by  magnetostatic  fields  will  have  profound  effects  on  resonance  phenomena  in  thin 

51 

films.  Such  effects  have  been  observed  by  Comly,  Penney,  and  Jones  at  high  microwave 
power  levels;  their  results  show  that  the  bottom  of  the  manifold  is  indeed  correctly  described 
by  Eq.(V-34a)  with  ^  =  0^^ 

Spin-wave  solutions  of  the  linearized  equations  of  motion  (V-25a,b)  and  (V-29a,b)  also  exist 

about  unstable  equilibria  of  the  uniform  mode  {4^^,  which  are  associated  with  magnetization 

reversal  processes,  and  which  may  be  found  from  Eqs.  (V-32a-c)  with  the  inequality  of  (V-32c) 

reversed.  The  dynamic  trajectory  of  the  uniform  mode  may  not  pass  through  the  point  4^^  =  ^ » 

but  because  of  the  condition  (V-26)  [see  (IV-51)  and  (IV-24)]  this  makes  very  little  difference  in 

the  dispersion  relation  for  a;_^.  It  has  been  pointed  out^^*^^"^"^  that  the  spin-wave  eigenfrequen- 

k 

cies  will  then  be  imaginary  for  a  certain  set  of  wave  vectors,  and  that  unstable  solutions  with 
time  dependence 
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Fig.  8.  Resonant  modes  af  a  thin  film  in  a  planar  field  of  800  oe  in  thin-film 
approximation  (average  aver  thickness),  spin-wave  approximation  (neglect 
magnetostatic  boundory  conditions),  and  mognetostatic  approximation  (neg¬ 
lect  exchange  fields)  (for  a  TPF). 
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Fig.  9.  Schematic  illustration  of  inverted 
and  normal  planar  mode  spectra  (thin-film 
approximation). 
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e 


(V-36) 


n_^t 

^  (n  =  iw  >  0) 
k  k 

may  be  important  in  switching.  We  now  show  that  the  growth  time  (for  increase  in  amplitude  by 
a  factor  e)  of  the  fastest  growing  spin  waves  is  approximately  the  same  as  the  coherent  rotation 
reversal  time,  so  that  switching  is  completed  before  any  spin-wave  amplitudes  increase  much 
above  their  equilibrium  values. 

It  is  easy  to  see  from  the  dispersion  relation  {V-33a)  (remembering  that  we  are  only  inter- 

2 

ested  in  low  fields,  i.e.,  cj,  «  cj  ,  cj  «  cj  )  that  unstable  spin  waves  (cj  <  0)  will  exist 

'  h  m  a  m'  ^  k 

providing 

2  A  A  A 

-(O)^  sin  +  “  J^X)  <  cos{cp^  -  /3)  +  cos  Zcp^ 

+  0)  L^k^<-a)  xsin^(^  -  ^)  .  (V-37) 

e  m'^  '^o  '  ' 

The  fastest  growing  of  these  will  have  ^  (propagation  along  m^)  and  kL  «  1,  so  that  their 

growth  rates  are 

fi  w  r-o),  cos($  -/3)-a)  cos  2^  -  w  >  0  (V-38) 

m  ^  h  '^o  a  ^o  e  ^  ' 

In  Fig.  9  we  sketch  the  ’* inverted”  spectrum  (m^  unstable)  and,  for  comparison,  the  normal  spec¬ 
trum.  From  (V-38)  we  have  the  important  result  that  the  growth  rate  of  spin  waves  in  a  thin 
film  in  unstable  equilibrium  in  low  fields  («47rM^)  is  bounded  by  the  growth  rate  of  the  uniform 
mode,  or 


=  y(47rM^)^/^  [-Hp  cos(^5^  -/3p)  cos  cos  2$^]^/^  (V-39) 

(k=^0) 


where  we  have  generalized  the  external  field  H  as  described  on  p.34.  This  result,  which  strongly 
suggests  that  spin-wave  instabilities  will  not  greatly  influence  the  behavior  of  the  uniform  mode, 
may  be  made  more  precise  by  comparing  the  maximum  growth  rate  for  spin  waves,  11^,  with 
the  inverse  switching  time,  t^  obtained  in  Sec.  IV.  From  (IV-47),  and  using  the  approximation 
suggested  on  p.  42,  we  find 

=  y  [4-M^(Hp  -  (V-40) 


where  J  ~  1,  except  for  =  0  and  =  7r,t  and  reduces  to  l/[K(sin  2  )3)  ]  for  =  0.  The 

pulse  threshold  field  is  the  field  at  the  transition  from  stable  to  unstable  equilibrium,  or 


H 


.  cos  {(p  —  )3  )  +  H,  cos  (3  —  /?,  )  +  cos  23  =0 
pt  '^o  ^p'  b  ^o  ^b  K  ^o 


(V-41) 


Inserting  (V-41)  in  (V-39)  we  obtain 


t  P  =  IT  is  a  pathalogica!  case  of  balanced  unstable  equilibrium,  and  we  shall  exclude  it  from  aur  treatment. 
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(V-42) 


Q  =y(-4;rM(H  -  H  J  cos  ($  - 

O  ‘  o  p  pt  o  ^  ^ 

- 1 

which  is  never  much  greater  than  t^  . 

The  results  obtained  in  Sec.  V-C  are  valid  whatever  the  origin  of  the  magnetic  fluctuations 
—  thermal  agitation,  inhomogeneities,  anisotropy  dispersion,  etc.,  —  and  are,  therefore, 
completely  independent  of  our  random  anisotropy  model.  In  the  following  pages  we  shall  consider 
transient  spin  waves  in  a  step  function  field  taking  as  initial  values  the  dispersion-induced 
ripple  components  of  Sec.  111. 

D.  TRANSIENT  RESPONSE 

W'e  first  rewrite  (V-25a,  b)  with  the  normalized  variables  [see  {IV-22)  to  (lV-24)] 

T  =  /o),  cj  t  =  y  MttM  H  t  {V-43a) 

V  h  m  V  o  p 

(V-43b) 

where 


€  = 


=  0(10"^) 


{V-43c) 


These  variables  have  been  chosen  such  that  during  magnetization  reversal 

0(1),  as  may  be  inferred  from  (IV-38a,  b)  and  {IV-50).  {The  dot  indicates  differentiation  with 
respect  to  r.)  The  equations  of  motion  of  the  spin-wave  components  are  now 

e  =  -G  -  G  (p  (V-44a) 

k  k  k 


+  G 


0  (p 


k 


{V-44b) 


with  initial  conditions 

e_(0)  =  0  {V-45a) 

k 

given  by  (111-35)  {V-45b) 

k 


where 

^ee  “  ^  'x  +  e  {CTq  [x  cos^  ((;()^-4>)  -  x]+  (AXg  k^  +  sin^  (p ^  ^  -  +  0(  e^) 

^  sin^(<p^ -4>)  +  e  [  AX^  k^h^  ^ -p '((p^)  ]  +  0(  e^) 

G_  =G  = -y  a  X  sin  2(<p  -<I>)  +  CXe^) 

Q(p  (pQ  2  o  o 


(V-46a) 

(V-46b) 

(V-46c) 


Here  we  have  written  the  exchange  field  in  terms  of  the  exchange  length  X^,  defined  by  (111-36). 
The  remaining  parameters  in  (V-46a,  b)  are 

A(^b*  2<;p^(0)  +  h^  cos  [/3^  -  <;^>q(0)  I  (V-47a) 
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(i  =  p,  b,  etc.) 


(V-47b) 


H. 


h.  = 


i  ~  H 


K 


b  1 

-P'(<Pq)  =  cos(<p^  “  IT  ~  h~ 

^  P  P 


(V-47c) 


(see  (III-32)  and  (IV-21c)].  The  uniform  mode  is  described  by  (IV-27a,b)  without  damping,  or 


(p  -  (j 

o 


(V-48a) 


•  b  1 

(j  =  p(<27  )  =  — sin((27  —B  )  — —  sin((p  —  J3,  )  —  sin2(p 

o  ^o  ^o  '^p  h  ^o  ^b  2h  ^o 

^  P  P 


with  initial  conditions 


(V-48b) 


h^  sin[<p^(0)  —  {^^]  +  2  sin2(p^(0)  =  0 


CT^(O)  =  0 


(V-48c) 

(V-48d) 


so  that 


p(<Pq(0)1  =  sin  |/3p  -  <p^(0)]  s  Pq  >  0 


For  the  case  (IV-19a,  c)  of  a  hard  axis  bias  field  and  an  easy  axis  pulse  field, 

1 

p{<p^)  =  sin  +  —  cos  v’o  -  ^  sin  Z<p^ 

Po  =  ^i  ■ 


(V-48e) 

(V-49a) 

(V-49b) 


The  finite  rise  time  of  the  pulse  field  will  not  be  considered  in  this  treatment;  thus,  p(t)  =  0  for 

T  <  0. 

We  next  eliminate  ©_^  from  (V-44a,b),  obtaining 
k 


00  .  ^ 
-  r — 
k  ^00  k 


-  G^  +  G 


00 


00  (p(p  Q(p  Q(p  G 


~  G 


00 


Q(p 


(p^  -  0 

k 


with  initial  conditions 


ipJO)  =  0 
k 

k 


(V-50a) 

(V-50b) 

(V-50c) 


As  can  be  seen  from  (III-35),  and  as  discussed  in  Sec.  III-D,  short  wavelength  ripple  components 
are  suppressed  by  exchange  forces.  We  therefore  confine  ourselves  to  the  long  wavelength  end 
of  the  spectrum 


k  <  X 


-1 


(V-51) 


and  in  Sec.  VI  we  will  show  that  the  main  contribution  to  the  spin-wave  torque  on  the  uniform 
mode  comes  from  components  with  k  ~  With  this  restriction,  the  coefficients  of  powers 

of  €  in  the  of  (V-46a-c)  (and  their  time  derivatives)  are  all  CXl). 
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We  are  now  in  a  position  to  greatly  simplify  Eq.  (V-50a).  Inserting  the  coefficients  G.^ 
from  (V-46a-c),  and  dropping  negligible  terms,  we  find 

<p^  +  €^x~^a  (r)  +  €”^xx  sin^  -  4>  ]  =  0  (V-52) 

k  k  ^  k 

where  a^(T)  =  0(1).  We  next  define  a  large  parameter  X  by 


X  = 


(V-53) 


For  kL  «  1 

X  =  (1  ~  I  kL  +  .  .  .) 


(V-53-1) 


and  for  a  TPF  with  k  =  A  ^  and  H  =5  oe,  A  20. 

e  p 

(p  +  A  ^  a^(T)  ip  +  (Af(T)]^  (p  =  0 
k  ^  k  k 


Equation  (V“52)  may  now  be  written 


(V-54) 


where 


a^fr)  =  xa^(T)  (ja^l-^laj)  (V-55a) 

f(T)  =  sin  [<;f>^(T)  -  4> )  .  (V-55b) 


-2 

Clearly  the  "damping”  term  A  a^<^  is  negligible  compared  to  the  "potential” 

k 

since  \ip  |  ~  A  |  f</?  |.  With  the  notation 

k  k 


term  A 


2^2 


k 


<p^{r) 


=  ^^(0) 

(V-56) 

k 

our  problem  is  now  reduced  to  solving  the  equation 

ip  +  [Af(T)]^  (p  =  0 

(V-57) 

with  initial  conditions 

II 

o 

(V-58a) 

<p(0)  =  0 

(V-58b) 

where  f(T)  is  given  by  (V-55b),  (p  is  determined  by  (V-48a-e),  and  A  »  1. 

2  ^ 

The  potential  (Af)  has  arisen  essentially  from  volume  poles  which,  in  equilibrium,  strongly 
attenuate  ripple  components  propagating  in  all  directions  except  along  we  recall  from 

Sec.  III-D  that  the  zone  of  nonattenuated  components  is  extremely  narrow.  (See  Fig.  5.)  In  fast 
magnetization  reversal,  however,  volume  magnetostatic  fields  quickly  build  up  to  large  values, 
since  components  initially  propagating  along  rn^  ”have  the  rug  pulled  out  from  under  them” 
when  m^  rotates,  and  go  through  a  transient  state  resembling  the  energetically  unfavorable 
transverse  ripple.  This  is  shown  schematically  in  Fig.  10.  On  the  other  hand,  although  the 
initial  very  small  amplitude  components  with  k  in  all  directions  other  than  =*=  do  go  through 
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INITIAL  TRANSIENT 

Fig.  10,  Schemotic  illustrotion  of  initiol  ond  transient  spin-wove  stotes 
(LMR)  in  rotationol  switching,  showing  creotion  of  volume  poles. 

a  transient  state  of  low  magnetostatic  energy,  their  amplitudes  will  not  grow  much  above 

equilibrium  values  unless  very  fast  relaxation  can  occur  (relaxation  time  «  t  ).  However, 

there  is  no  evidence  for  such  a  rapid  process,  and  in  any  event  we  are  neglecting  all  relaxation 

in  this  treatment.  Thus,  we  can  conclude  that  the  primary  effects  in  magnetization  reversal  will 

come  from  spin  waves  with  4>  <2?  (0)  and  ^  ^  cp  (0)  tt,  that  is,  the  initial  LMR.  But  since 

^  ^o  o 

6M  is  real,  we  see  from  (V-5a)  that 

m  _  =  m *  (V-59) 

-k  k 

and  it  is  therefore  sufficient  to  solve  (V-57)  for  4> 

The  presence  of  the  large  parameter  \  in  (V-57),  coupled  with  the  fact  that  f(T)  is  not 
rapidly  varying,  suggests  the  WKB  method,  for  which  we  refer  to  Appendix  C,  Our  main  interest 
is  in  the  asymptotic  solution  (C-14),  which  is  valid  except  near  free  points,  at  which  Rt)  =  0, 

The  magnitude  and  phase  of  this  solution  is  determined  from  the  initial  conditions  (V-58a,b); 
but  for  the  components  of  interest,  r  =  0  is  at  or  very  near  a  free  point.  Therefore,  a  solution 
joining  the  asymptotic  region  to  the  neighborhood  of  the  free  point  must  be  found.  W’e  first 
make  the  substitutions 


t](t)  =  <pAt)  -  (P  (0) 

0  0 

(V-60a) 

=  4>  —  cp  (0) 

0 

(V-60b) 

so  that  (V-55b)  becomes 

Rr)  =  sin  (77(7)  -  4'] 

(V-61) 

with  t]{t)  determined  from  (V-48a-e),  which  become 

7?(t)  = 

(V-62a) 

o‘q(t)  =  p(t]) 

(V-62b) 

T7(0)  =  0 

(V-62c) 

(T  (0)  =  0 

0 

(V-62d) 

p(0)  =  P^  >  0  . 

(V-62e) 

The  asymptotic  generating  function  for  [Xf(T)]^  is  given  by  (C-15) 
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no  d? 


(V-63) 


8(0  .  X  f 


where  is  a  constant,  as  yet  unspecified.  We  examine  g  in  the  vicinity  of  the  free  point  77  =  'I', 
with 

I'l'  I  «  1  (V-6‘1) 

in  accord  with  the  previous  discussion. 

Consider  first  the  case  'I'  ^0,  for  which  we  let 


(V-65) 


Since  j]  is  initially  zero  and  xt  is  small,  the  free  point  may  be  found  by  expanding  r](T)  in  a 
Taylor's  series  about  r  =  0,  with  the  result 


provided 


'I'  « 


Now  we  let 


X 


(V-66a) 


(V-66b) 


(V-67) 


and  expand  f  about  obtaining 


where 


.2  3  4 

f(x)  -  p  It  X  +  ~  +  P.,T  +  p.  tt  ^ 
o  \  o  2!  '^1  o  31  ^14! 


P^  =  P’(0) 


) 


(V-68a) 


(V-68b) 


Inserting  this  result  into  (V-63),  we  find 

_Z  3 


g(x)  =  xp^(t^  it  ^  fr  +  '^ro  It  ^  '’i  fr  ^ 


(V-69) 


For  the  limiting  case  xt  =  0  (t  =  0),  the  x  term  of  (V-69)  will  approximate  g  for  x  in  the  range 
S  given  by 


x^  « 


(V-70) 


Then,  following  the  argument  of  pp.79-8  0,  we  find  that  the  approximate  solution  to(V-57)  in  is 

The  asymptotic  expansion  of  <p  [given  by  (C-11)  and  valid  in  defined  by  (C-lOb)]  is  valid  in 
part  of  inequality  (C-23)  is  satisfied  for  some  x^  in  with  r  =  l/6  and  f^  =  P^/^, 

(C-23)  becomes 
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(V-72) 


Then,  combining  (V-70)  and  (V-72)  we  arrive  at  the  condition  for  overlap  of  ^oo 


X  » 


(^r 


(V-73) 


which  is  well  satisfied  except  as  the  excluded  case  of  a  pulse  field  antiparallel  to  m^(0). 

(See  footnote  on  p.  51.) 

Now,  suppose  ^  >  0.  Under  what  conditions  will  the  solution  (V-71)  still  join  the  asymptotic 
region  to  the  vicinity  of  the  initial  point,  x  =  First,  (V-73)  must  be  satisfied;  but, 

second,  and  more  important,  the  largest  value  of  |x|  for  which  the  x^  term  of  (V-69)  is  not 
negligible  compared  to  the  x  term  must  be  less  than  the  smallest  value  of  |x|  in  .  F'or, 
otherwise,  the  asymptotic  solution  will  join  not  with  (V-71)  but  with  a  solution  involving  also 
(or  exclusively)  the  first  term  of  (V-68a).  This  condition  is  met  if 

(V-74) 


as  may  be  seen  from  (V-69)  and  (V-72);  using  (V-66a),  (V-74)  becomes 


«  yp  =4-  «  1 

2  '^o  o 


(V-74-1) 


If  4^  <  0,  as  given  by  (V-66a)  is  imaginary;  there  is  a  nearly  free  point  at  r  =0,  but  no 
free  point  near  77(7)  =  0  for  real  r.  However,  we  may  consider  x,  defined  by  (V-67),  as  a 
complex  variable  and  by  analytic  continuation  extend  the  definitions  of  f(x)  and  g(x)  into  the 
complex  domain.  Then  all  the  results  for  4^  >  0  are  equally  valid  for  4'  <  0,  provided  x,  r^, 
and  ^  are  replaced  by  their  absolute  values  in  (V-66b),  (V-70),  (V-74),  and  (V-74-l). 

Collecting  our  results,  we  have  found  that  a  solution  to  (V-57),  valid  from  the  neighborhood 
of  T  =  0  up  to,  but  not  including,  the  next  free  point  (near  17  =  tt),  is 


provided  the  angle  4^  between  the  initial  direction  of  m^  and  the  wave  vector  k  of  the  component 

(p  satisfies 
k 

|4^  I  «  4^^  «  1  .  (V-76) 


Note  that  what  we  have  found  is  essentially  a  4^  =  0  solution,  as  it  is  based  on  the  4^  =  0  generating 

function  [(V-69)  with  =  0].  For  14^1  »  4^^,  it  is  not  hard  to  show  that  a  solution  of  the  form 

(V-75)  exists  with  r  =  ^.  But  for  |  4^  |  ~  4^  ,  no  such  solution  is  possible  because  the  asymptotic 

°  2  3 

form  must  be  connected  to  a  solution  in  a  range  of  x  for  which  the  x  and  x  terms  of  (V-69) 
are  the  same  order  of  magnitude.  However,  we  shall  show  in  Sec.  VI  that  the  strong  inequality 
(V-76)  does  not  exclude  any  physically  important  components  of  the  ripple. 

The  constants  are  evaluated  from  the  initial  conditions  (V-58a,b).  At  r  =  0,  x  =  — r^, 
so  that  from  (V-69)  and  (V-74) 


lglT=0 


(V-77) 
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Then,  using  the  expansion  of  J^(g)  about  g  =  0/ 
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Lj  ^2^  m!  r(y  +  m  +  1) 
m=0 


and  retaining  only  the  leading  term  of  g 

Ap  - 
/  \  o  3 

g(x)  X 


(V-78) 


(V-79) 


we  find 


<p(x) 


yi 


■  B 

-1/6  B^ 

1 

lo 

h 

r(|) 

\  12 ) 

\\Z  )  ’ 

(V-80) 


Then,  since  at  r  =  0 


|x|  «  (\p^) 


-1/3 


(V-81) 


the  initial  conditions  (V-58a,b)  lead  to  the  result 


B_  -  r(|)  (|)^/^ 


(V-82a) 

(V-82b) 


Finally,  we  obtain  the  asymptotic  expansion  of  ip(r) .  The  asymptotic  region  is  given  by 
(C-lOb) 


or  [see  (V-79)  ] 


(|  «T  <T^ 


(V-83a) 


(V-83b) 


where 


tt  +  ^ 


(V-83c) 


The  upper  bound  =  O(t^),  and  since  we  will  not  be  investigating  the  behavior  of  spin  waves  at 
the  end  of  the  reversal  process,  there  is  no  need  to  specify  r  ^  precisely.  From  (C-11),  noting 
that  </>  =  0  in  [R  ,  we  have  the  resultl^ 


T  in  S  ‘ 


\  ^  1 
T-)  -  xl 


(V-84) 


t  Except  for  chonges  in  the  constont  B  ond  the  phose  11/6,  (V-84)  Is  the  osymptotic  volue  of  </7(t)  for  oil  4^.  A 
computer  solution  of  (V-57),  however,  hos  shown  thot  the  true  constont  is  opproximotely  equol  to  B  for  , 

with  moximum  error,  ot  of  15  percent,  ° 
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Fig.  11.  Transient  response  of  a  planar  spin-wave  component  in  rotational  switching, 
or  [using  ( V-56)  ] 


T  in  31  ’  [f(T))‘^/^  cos  [g(T)  -  -g-] 

k  oo  k 


where 


Jl-r'(f) 


(V-85) 


(V-86) 


Figure  11  is  a  sketch  of  <p{t) .  The  spin-wave  component  normal  to  the  film  plane  may  be  found 
from  (V-44b)  and  (V-56): 


e_,  =  g~q  [<p(t)  -  Gp,^^(p(T)]  (p^(o) 


Q(p^ 


(V-87) 


and  with  the  help  of  (V-46a,  c)  and  (V-53)  we  find 


(V-88) 


Equations  (V-85)  and  (V-88)  comprise  the  main  result  of  Sec.  V-D.  They  show  that  except 
at  the  beginning  and  end  of  the  reversal  process,  =  0  spin  waves  (LMR)  precess  in  an  ellipti¬ 
cal  path  about  rn^(T)  at  a  frequency 


a> 

k 


g(T) 

T 


sinTj(0  d^ 


The  ellipticity  is  given  by 


(V-89) 
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^^{t)  = 
k 


^  sin7](T) 


(V-90) 


—  cot(g  -  t) 


- 1  /  9-1 

[In  a  TPF  the  important  set  of  components  with  k  ~  A  have,  at  77(7)  =  7r/2,  cj  ~  35  X  10  see 

°  k 

and  ~  |.]  We  note  that  the  results  we  have  obtained  depend  only  implicitly  on  the  uniform 
k 

reversal  mode,  except  at  the  very  start  of  switching.  Thus  even  though  we  have  neglected  damp¬ 
ing,  any  real  loss  mechanism  will  require  a  finite  time  to  become  effective,  and  therefore  may 
be  included  through  its  effect  on  particular,  the  perturbation  of  the  uniform  mode  by 

spin  waves,  which  will  be  treated  in  the  next  section,  need  not  be  small  for  the  theory  to  be 
valid;  be  determined  self-consistently  from  (V-6a),  (V-85),  and  (V-88). 
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VI.  EFFECT  OF  SPIN  WAVES  ON  MAGNETIZATION  REVERSAL 


A.  INTRODUCTION 

Although  in  static  equilibrium  ripple  has  very  little  influence  on  the  mean  magnetization 
m^,  the  dynamic  situation  is  quite  different.  As  noted  in  Sec.  V,  volume  magnetostatic  fields 
attain  large  values  during  fast  rotational  switching  because  k  (for  the  largest  amplitude  compo¬ 
nents)  and  m^  no  longer  remain  parallel.  From  an  energy  viewpoint,  this  transient  transverse- 
component  ripple  is  created  at  the  expense  of  the  energy •  m^(0)  which  the  uniform  mode  has 
acquired  from  the  pulse  field.  The  transient  spin-wave  state  therefore  provides  a  sought-after 
dissipation  mechanism  for  m^,  at  least  during  the  first  quadrant  of  switching  [0  t](t)  ^  7r/2, 

where  17(7)  =  <P When  t}  >  tt/Z,  if  these  components  cannot  lose  their  energy  to  the 
lattice  or  to  other  components  they  will  return  it  to  the  uniform  mode,  thereby  conserving  energy 
(since  magnetostatic  fields  are  decreasing).  Thus  we  expect  damping  in  the  first  quadrant  and 
anti-damping  in  the  second;  this  is  indeed  what  we  shall  find  in  the  calculation  to  follow. 

However,  a  much  more  striking  effect  is  anticipated;  spin-wave  locking  of  the  uniform  re¬ 
versal  mode,  which  will  occur  if  the  magnetostatic  energy  of  the  intermediate  state  is  greater 
than  the  energy  available  for  rotational  switching.  Since  this  latter  energy  decreases  with  de¬ 
creasing  Hp,  we  may  expect  to  find  a  critical  field  below  which  the  uniform  mode  is  locked. 

In  Sec.  VI-C,  we  calculate  this  critical  field,  which  provides  an  important  point  of  contact  with 

experiment.  We  also  suggest  how  rotational  switching  might  take  place  with  H  <  H 

p  pc 

These  phenomena  r'ny  alno  unHerstood  from  a  torque  viewpoint.  The  nonlinear  spin-wave 
reaction  torque  on  m^,  which  we  shall  find  to  be  predominantly  z-directed,  retards  reversal  in 
the  first  quadrant  and  accelerates  it  in  the  second.  If  at  any  point  in  the  reversal  process  with 
m^  in  the  first  quadrant  the  reaction  torque  becomes  equal  in  absolute  magnitude  to  the  uniform 
reversing  torque  i^  •  m^  x  h^,  the  uniform  reversal  mode  will  become  locked.  We  now  put  these 
ideas  on  a  quantitative  basis  by  computing  the  reaction  torque,  using  the  results  of  Sees.  Ill  and  V. 

B.  NONLINEAR  REACTION  TORQUE 

The  spatially  varying  part  of  the  magnetization,  6M(r,  t),  will  react  on  the  mean  magnetiza¬ 
tion  through  the  spatial  average  of  the  ensemble  average  of  the  nonlinear  contribution  to  the 
effective  torque.  (By  nonlinear,  we  mean  second  order  in  the  small  quantities  6M/M^  and 
6H/M^;  there  is,  of  course,  no  linear  reaction  torque.)  This  appears  in  the  equation  of  motion 
for  ^0^(1)  [see  (V-4a)]  as  a  reaction  torque 


(VI-1) 


or,  in  terms  of  the  Fourier  components  of  6M  and  6H  [see  (V-6a)], 


(VI-1-1) 


k^O 


We  recall  that  the  bar  denotes  an  ensemble  average  and  the  bracket  a  spatial  average.  The 
spatial  components  of  m_^  are 


k 


(VI-2) 


m  =M(i0  +10) 

k  °  ^  k  ^  k 

while  the  spatial  components  of  h  ^  are  found  through  (V-7)  as  linear  funetions  of  0  ^  and 

“k  -k 

(p  Using  (V-59),  we  see  that  m_^  ^  h  will  contain  terms  of  the  following  four  types; 

-k  k  -k 


■  (vi-3) 

k  k  k  k  k  k 

For  T  in  !K  (to  whieh  we  shall  confine  ourselves  in  this  section),  we  have  seen  that  the 

OO 

spin-wave  eomponents  (P_^{t)  and  0^(t)  oseillate  at  a  frequency  a;^(T)  much  greater  than  the 

k  k  ^  k 

precessional  frequency  of  m^.  Therefore,  m^  will  respond,  not  to  the  instantaneous  reaetion 

torque,  but  to  an  effective  torque  whieh  is  averaged  over  the  period  of  this  oscillation.  Since 

f(T)  is  slowly  varying,  we  may  consider  it  to  be  constant  during  one  period  of  oscillation.  Then 

since  and  (p^  are  90°  out  of  phase,  terms  of  the  last  two  types  in  (VI-3)  give  no  contribution 
k  k 

to  the  effeetive  torque,  and  we  are  left  with 

m  Xh  =M  It  K®‘^(|0  \^  -  \  0  |^)+T_K'^“’|<p  |^-Tk'^®|0  1^1  (VI-4) 

^  ol  ^  k  k  k  ^  k  k  ^  k  k  J 

2  2  2  I  2 

where  the  factors  sin  (g  —  7r/6)  and  eos  (g  — 7r/6)  of  |0^|  and  |<^^|  >  respectively,  are  to 

k  k 

be  replaeed  by  their  average  value  \  .  We  note  that  exchange  fields  (whieh  appear  in  the  diagonal 

elements  of  K  )  do  not  eontribute  to  the  effective  reaetion  torque.  Also,  since  «  47rM  V 
k  K  o 

[except  for  extremely  long  wavelength  components  which,  we  will  find,  make  no  significant 

contribution  (t  )  ],  the  anisotropy  term  in  K  ^  niay  be  safely  neglected,  leaving  only  magneto- 
^  ^  k 

static  terms.  Equations  (VI-4),  (V-22),  (V-23a,b),  and  (V-26)  then  give 

X  h  =  47rM^  7  ^(v  -  -!•)  (  ^  -  I  6^1  ^) 

k  -k  o  m  o  ^  o  ^  ^ 

+  ‘i0X  i  sin2(<p^-4.)  ^  [1  -X  sin^(<p^-4.)l  |0_^|^)  (VI-5) 

k  ^  k 

We  now  show  that  only  the  i^-eomponent  of  have  a  physically  important  effect  on 

m  .  Taking  the  ensemble  average  of  m  x  h  and  performing  the  sum  over  k,  we  write 
°  k  -k 


( t  )  =  47rM  {  i  lb  s  +  i^s^  +  i  ib  s  ) 

o  r  o  m^o  m  0  0  (p^o  (p 


where^ 


(VI-6) 


=  7  Tj  x(kU  sin2(«,^-4.)(|<^_|2-  |0_|2) 

'  k  k  ' 

k^O 


(VI-7a) 


t  Strictly  specking,  the  ensemble  averages  in  (VI-7)  should  Include  the  factors  containing  (p^.  However,  as  was 
noted  on  p.  21 ,  rondom  fluctuotions  of  (p^  vonish  in  the  limit  S  ^  ond  are  therefore  of  no  physical  importonce. 
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■  i  2  x(kL)  Sin2(<p^  -  #) 

iT^o 


(VI-7b) 


^  (1  -  X  (kL)  sin^  ((p  -  <!.)]  |ej^ 

^  k 

k^O 

From  (I-la),  (V-5a),  and  (VI-2)  we  obtain  the  condition 


<|6M|^> 

M  ^ 
o 


from  which  we  deduce  the  separate  conditions 


E  w^|ej2«i 
k  k 

k^O 


(VI-7c) 


(VI-8) 


(VI-9a) 


E  w  \(p  \^  «  i  (VI-9b) 

k  k 

k?^0 


where  w  is  any  weighting  function  satisfying 
k 


|w  I  <  1  (VI-9c) 

k 

[Note  that  with  w_^  =  1,  (VI-9b)  leads  to  the  restriction  (III-58)  on  the  magnetization  dispersion.) 
^  k 

Since  0  <x  4  these  conditions  may  be  applied  to  the  sums  of  (VI-7),  and  we  see  that 

|Sj(<Po)l  «  ^  j  =  m,  e,  .  (VI-10) 

Next  we  insert  the  effective  reaction  torque  of  (VI-6)  into  the  equations  of  motion  (V-4a)  for 

m  and,  with  the  help  of  (V-43a-c)  and  (V-48a,  b)  obtain 
o 

m 


_ o 

m 

o 


=  —a  s  ((p  ) 
o  m  ^o 


(Vl-lla) 


CT  =  p(<pj  -  e"^s-(<p  )  (Vl-llb) 

O  O  t7  o 

(b  =  a  [1  -  s  ((^  )]  (VI-llc) 

where  we  have  neglected  random  fluctuations  of  m^  and  a^,  as  well  as  cp (See  footnote,  p.  62.) 
We  note  first,  from  (VI- 11a),  that  |rn^|  =  longer  a  constant  of  the  motion  (|m|  = 

still  is,  but  now  and  m^  differ  by  (|6M|^)).  However,  the  relative  change  in  m^  during 
switching  is  small  (since  I  |  «  1),  and  no  profound  effects  are  anticipated.  We  next  note, 
from  (VI-llc),  that  the  direct  effect  of  spin  waves  on  (p ^  is  small  (and  therefore  of  no  physical 
importance);  with  the  help  of  (VI-llc)  we  may  integrate  (VI- 11a),  obtaining  the  implicit  solution 
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(VI-12) 


m  (t)  =  m  (0)  exp 
o  o 


(T) 

s  (4)d? 

(0) 


However,  because  the  change  in  is  probably  not  measurable,  we  shall  carry  this  solution  no 

further  (although  evaluation  of  s^  is  straightforward)  and  turn  instead  to  the  one  potentially  large 

spin-wave  effect  —  that  of  s^, 

-2  4  ^ 

Since  €  »  1  ('^10  ),  the  second  term  on  the  right  side  of  (Vl-llb),  which  comes  from  the 

i^-component  of  the  reaction  torque,  is  not  a  priori  negligible  and,  in  fact,  could  even  dominate 
the  uniform  field  torque  term  picp .  To  compute  s^  we  take  as  given  by  (V-85)  and,  re¬ 
calling  the  definition  of  X  (V-53),  obtain  ^ 


"  2  Po'> 


1,2  ,1/3^2  Y  ~2/3  -1/3 


C“  E  X  ■  X 
k^O 


COS  {(p  —  ^)  \(p  (0)  I 

°  k 


(VI-13) 


The  initial  ripple  amplitudes  are  given  by  (III-39)  and,  converting  the  sum  to  an  integral  [see 
(III-44)  I,  we  find 


,1/3  ^  1  C°°  ~Z/i  -1/3,,  ,  2,  2, -3/2 

o^  ^  j  Ml  +  r^k  ) 


Se(<Po>  =  p 


X  kdk 


(p  {0)+7r/2 


Mo 

/( 


'cpA0)-n/Z 


COS  {(p  —  4>)  d4) 


(VI-14) 


We  have  used  (V-59)  to  relate  (p  ^  to  cp^,  and  it  must  be  remembered  that  the  solution  for  <p^{r) 

-k  k  k 

we  are  integrating  is  valid  only  for  |  4^  |  <<  4^^.  The  second  integral  in  (VI- 14)  is 


tt/Z 


( 1  +  \  ^k^  +  \  L  sin^  4^)  ^  cos  (t](t)  —  4^  ]  d4^  i  "il'? 

7r/2  ^  ^  ^  (a-^b)^/ 


(VI-15) 


where 


2  2 

a  =  1  +  X  k 
e 


(VI-l6a) 


b  ~  X  I-j  ^  V 
m  ^ 


(VI- 16b) 


and  we  have  made  use  of  the  fact  that  for  k  =  0(X  ^),  b  »  a.  The  main  contribution  to  I,  is 

?  e  *  4#  . 

from  4^  b/a  <  1;  using  (III-37),  (V-74-1),  and  (V-53),  we  see  that  for  components  with  k  =  ) 


2b  p  .  ,2/3 

'^o  ¥  8]^  <^Po) 


(VI-17) 


-1 


Since  x.  A,  and  h^  are  all  ~1,  if  X  »  (we  recall  that  X  »  1  and  the  limit  0  is  excluded), 

then 


o  a 


(VI-18) 


and  only  components  for  which  our  solution  is  valid  contribute  to  I,  .  But  for  typical  values  of 

2 

the  parameters  involved,  (VI- 18)  is  satisfied  only  marginally,  e.g.,  4'^b/a  ^  2.  However,  as 
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noted  in  the  footnote  on  p.  58,  the  solution  may  be  extended  to  |^  |  <  without  serious  error. 
We  estimate  the  possible  error  in  to  be  +20  percent. 

Equation  (VI- 14)  now  becomes 


7  Pr 


(VI-19a) 


where 


,  r"  ~l/6  -1/3  ,,,  ^  2,2.  ,,  2.2.  ,-3/2 

1^=3  X  X  ((1  +  )  (1  +  Xgk  )  ] 


kdk 


(VI- 19b) 


We  evaluate  by  the  same  method  used  to  find  6  on  pp.  27-29,  with  the  result 


l,"/‘ 


dk 


(VI-20) 


where  I(p,  q)  is  given  by  (III-65).  Thus 

.  ,2„, 13 


„  _  ,  2  1/3  r(  12)  /rr\2  3/  2/3 

''  ■  '  V4A  )  ^ 


COS7/(t)  .  (VI-21) 


It  is  convenient  to  express  s^  in  terms  of  the  zero  field  magnetization  dispersion  6^,  given  by 
(III-66)  or  (III-67)  with  A  =  1,  and  we  finally  obtain  for  the  reaction  torque  contribution  to  d 


-2  .  -2/3,  „  ,  , 

c  s^  =  h  '  d  R  cos  t](t) 
Up  o 


o 

(VI-22a) 


where 


o  ^o 


■Uo*/  ^ 


R  =  C 


(VI-22b) 


(VI-22C) 


and 


c^  =  (|)'/'  r(||)  r(4) 


r<6> 


Lr(i) 


=  2.45 


(VI-22d) 


Before  we  investigate  the  effect  of  this  reaction  torque  on  the  uniform  mode,  it  is  helpful 

to  recall  the  meaning  of  the  factors  in  (VI-22a).  h^  is  the  (constant)  magnitude  of  the  pulse  field 

H  .  normalized  to  the  uniform  anisotropy  field  =  2K  /M  .  The  factor  d  is  a  function  of 
p*  K  ^  o'  o  o 

the  initial  state  of  the  film  and  of  the  direction  of  H  ,  and  with  and  A  given  explicitly  by 

(V-48e)  and  (III-32)  is 


=  {sin  (Pp  -  (p^(0)  (cos  Z(p ^{0)  +  h^^  cos  -  (p^(0)  ]}" 

For  the  case  of  a  hard-axis  bias  field  and  an  easy-axis  pulse  field  (see  (IV-19a,  c) 


X-5/12 


(VI-23) 
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(VI-23-1) 


d 

o 


( 1  - 


-5/12 


and  for  the  case  of  zero  bias  (see  (IV-l9b)] 

d  =  (sin 
o 


(VI-23-2) 


The  reaction  torque  coefficient  R  is  a  function  of  the  thickness  and  magnetic  parameters  of  the 
film;  for  a  TPF,  R  =  456^  =  0.26.  The  last  factor,  cos  17(7)  =  cos  —  </?^(0)  ],  introduces 

an  implicit  time  dependence  into  the  reaction  torque.  The  change  in  sign  of  the  torque  at  17  =  7/2 
is  expected,  since  for  k  ||  *^^(0)  tbe  magnetostatic  field  and  planar  spin-wave  components  are 
parallel,  producing  no  reaction  torque. 


C.  ROTATIONAL  MAGNETIZATION  REVERSAL:  DYNAMIC  LOCKING 

The  equations  of  motion  of  the  uniform  mode,  as  modified  by  the  ripple,  are  [from  (VI-1  lb,  c) 
and  (VI-22a)  ] 


<^q(t)  ^  p[t7(t)]  -  V  cost7(t) 
77(t)  =  a^(T) 


(VI-24a) 

(VI-24b) 


where  the  uniform  torque  [from  (V-48b)  ]  is 

h,  . 

p(t7)  =  -sin  [t7  4  ^  y-  sin  [t?  +  <P ^  sin2  [t?  +  (0)  ]  (VI-24c) 

P  P 

and 

V  =  R  (VI-24d) 

p  o 

providing  r  is  in  given  by  (V-83a-c).  (The  initial  reaction  torque  is  not  V,  but  zero,  since 
M  is  in  equilibrium.)  K  V  is  small  (’'small"  will  be  defined  below),  rather  than  solve  (VI-24a,  b) 
we  compare  them  to  the  corresponding  equations  with  phenomenological  damping,  (IV-27a,  b), 
and  verify  the  observation  made  in  Sec.  VI-A,  namely  that  the  reaction  torque  results  in  an 
effective  damping  in  the  first  quadrant  and  an  effective  anti-damping  in  the  second.  An  explicit 
solution  for  rfir)  would  then  show  retardation  in  the  first  quadrant  and  acceleration  in  the  second, 
with  a  net  effect  of  somewhat  slower  switching  since,  without  damping,  the  greater  portion  of 
the  time  for  reversal  is  spent  in  the  first  90®.  This  may  account  for  the  anomalous  damping 
observed  in  coherent  rotation  (see  p.  2);  however,  to  compare  theory  and  experiment  on  this 
point  the  intrinsic  damping  torque  should  be  included  on  the  right  side  of  (VI-24a),  and 

a  machine  computation  becomes  necessary. 

But  now  suppose  that 

V  cost7(t)  >  p|t7(t)1  for  0  <  t  ^  (VI-25) 

where  is  the  smallest  t  in  51^  (usually  «  i).  Then  the  effective  torque 

Peff(’?)  s  piv)  —  V  cos  t)  (VI-26) 
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becomes  negative  after  a  short  initial  buildup  time  for  the  reaction  torque,  and  rotational  reversal 
is  prevented;  the  uniform  mode  has  become  locked.  However,  the  inequality  (VI-25)  is  sufficient 
but  not  necessary  for  locking,  which  may  occur  later  in  the  reversal  process.  We  define  a  crit¬ 
ical  value  of  V,  V  ,  as  the  minimum  value  for  which  the  uniform  mode  becomes  locked;  this 
c  ' 

will  occur  at  an  angle  r} with  0  <  77^  <  tt/Z;  and  t} ^  may  be  found  most  easily  by  a  generali¬ 
zation  of  (IV -3 7): 

=  0  (VI-27a) 

V  =  0  (VI-27b) 


These  relations  determine  77^  and  either  V^  (at  fixed  field)  or,  alternatively,  the  critical  pulse 
field  hp^  for  irreversible  rotation  through  77^  (at  fixed  R) .  We  shall  adopt  the  second  viewpoint, 
since  the  field  is  conveniently  varied  experimentally.  With  the  help  of  (VI-24c)  and  (VI-26), 
(VI-27a,  b)  become 


1 

X  sin0  +  Y  cos  (p  —  sin2(z?  =  0 
^oc  ^  oc  2  ^  oc 

X  cos  <p  —  Y  sin  (p  —  cos  Z0  =0 
^  oc  ^  oc  oc 

where 

0  =77  +  (0) 

^oc  'c 

X  =  — h  cos/3  —  h,  cos /3,  —  h  d  R  sinc^  (0) 
pc  p  b  b  pc  o  ^o 

Y  =  h  sin/3  +  h,  sin/3,  —  h  d  R  cos  (p  (0) 
pc  ^p  b  ^b  pc  o  ^o 

We  find  from  (VI-28a,b)  the  relations 


(VI-28a) 

(VI-28b) 

(VI-29a) 

(VI-29b) 

(VI-29C) 


X  =  cos  0 

^  oc 


(VI-30a) 


Y  = 


.  3 

sin  0 
^  oc 


(VI-30b) 


provided  (p^^  >  from  which  we  obtain  the  equation  for  the  critical  field 


+  y2/3  1 


(VI-31) 


In  Appendix  D  we  examine  (VI-31)  for  two  field  configurations  of  experimental  interest 
((IV-19a,b)]. 

For  0  <  0  {t  ) .  locking  occurs  at  the  start  of  reversal  —  at  some  time  between  0  and  r  . 

Then  to  find  the  critical  field  we  would  compute  the  reaction  torque  from  spin  waves  given,  not 
asymptotically,  but  by  (V-75)  and  (V-82a,b).  We  shall  not  attempt  this  formidable  integration, 
but  instead  make  the  rather  crude  assumption  that  R  =  0  for  t  <  Then  in  the  limit  0, 

the  critical  field  is  determined  by  the  disappearance  of  the  initial  torque 

p^fj.(0)  =  0  (VI-32) 
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from  which  we  obtain 


h  =  d  R  CSC  W  -  (p  (0)  1 
pc  o  ^  p  ^o  ^ 


(VI-33) 


Note  that  (VI-33)  implies  a  cube-law  dependence  of  h^^  on  the  dispersion  6^,  which  should  be 
very  striking  if  it  actually  exists. 

We  conclude  with  a  general  discussion  of  the  critical  field  and  its  implications.  It  should 

first  be  noted  that  for  R  =  0,  h  becomes  h  ^  —  the  usual  threshold  field  for  irreversible  uniform 

pc  pt  ^  ^ 

rotation  [see  (V-41)],  and  (VI-31)  reduces  to  the  Stoner-Wohlfarth  asteroid 


(— h  .  cos  (3  —  h,  cos  ) 

pt  p  b  b' 


+  (h  ,  sin/3  +  h,  sin/3,  )^^^  =  1 
'  pt  '^p  b  '^b 


(VI-34) 


Then  from  (VI-29b,  c),  (VI-30a,  b),  and  (VI-31) 

ah 


pc 


aR 


1/3  cos  -  <p^(0)  ] 
R=0  pt  o  sin(^p-„^^) 


(VI-35) 


For  small  R,  therefore,  h^^  —  h^^  is  proportional  to  R,  and  from  (VI-22c) 


h  —  h  .  oc  (5 
pc  pt  o 


for  R  «  1 


(VI-36) 


We  next  note  that  for 


h  ^  h 
P  pc 


(VI-37) 


rotational  switching  is  modified  but  not  locked  by  the  ripple.  Thus  we  identify  the  field  region 
(VI-37)  with  that  of  high-speed  switching,  or  "coherent  rotation"  (region  3  of  FTg.  1),  and  the 
critical  field  h^^  with  the  threshold  field  for  coherent  rotation  (h^^  of  Fig.  1).  Preliminary 
experiments  performed  by  T.  D.  Rossing  and  G.  P.  Weiss  at  the  Laboratory  have  verified  the 
dependence  of  this  threshold  on  the  bias  field  h^,  as  given  by  (VI-31).  (See  Fig.  D-2.)  If  6^ 
can  be  measured  independently,  the  dependence  of  h^^  on  the  magnetization  dispersion,  as 
given  by  (VI-36)  or  its  equivalent  for  R  not  small  (see  Appendix  D),  should  provide  a  sensitive 
experimental  test  of  the  theory.  Analogous  with  the  replacement  suggested  on  p.  42  to  approxi¬ 
mately  account  for  anisotropy  effects,  we  may  replace  H,  in  the  solutions  found  in  Sec.  IV  for 
(p^  and  a^,  by  and  thereby  approximately  account  for  both  anisotropy  and  spin-wave 

effects . 

For 


pt 


<  h  <  h 
P  pc 


(VI-38) 


m^  starts  to  rotate  but  at  some  point  in  the  first  quadrant  becomes  dynamically  locked  by  the 
magnetostatic  field  of  spin  waves  propagating  in  the  direction  of  ^^^(0)  (LMR).  We  identify  the 
field  region  (VI-38)  with  that  of  intermediate-speed  switching,  or  "noncoherent  rotation," 

(region  2  of  Fig.  1),  and  offer  the  following  conjecture  on  the  reversal  process  following  locking: 
We  have  seen  that  the  occurrence  of  dynamic  locking  depended  on  the  relaxation  time  of  initial 
ripple  components  not  being  short  compared  to  the  switching  time.  However,  with  the  uniform 
mode  locked  these  components  can  now  relax,  either  to  the  lattice  or  to  components  propagating 
in  the  instantaneous  direction  of  m^,  thereby  losing  their  large  magnetostatic  energy  and  un¬ 
locking  the  uniform  mode.  Thus  we  envisage  a  highly  damped  rotational  process,  with  switching 
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speed  controlled  by  the  spin-wave  relaxation  rate,  and  with  a  threshold  field  h^^  ^ 

This  is  in  agreement  with  experimental  observations  (see  Sec.  I),  but  no  quantitative  assessment 
of  switching  speeds  expected  can  be  made  until  spin-wave  relaxation  processes  are  better  under¬ 
stood,  or  at  least  until  relaxation  times  have  been  measured. 
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APPENDIX  A 

NOTATION;  PARAMETERS  OF  A  TP F 


L  NOTATION 

The  following  subscripts  appear  on  H  (magnetic  field),  T  (torque  =  M  X  H),  and  their  Fourier 

components  h  and  t  : 

k  k 


Subscript 


Meaning 


eff  effective 

m  magnetostatic  (omitted  in  Sec.  11) 

e  exchange 

a  anisotropy 

h  external  field  (at  an  angle  (S  to  i^) 

p  external  pulse  field  (at  an  angle  to  i^) 


The  following  subscripts  appear  on  H  (external  field  magnitude),  fS  (angle  of  external  field 
to  i^),  and  H  external  field  normalized  to  uniform  uniaxial  anisotropy  field 

"k  = 


Subscript 


Meaning 


p  pulse 

b  bias 

11  pulse  along  — i^  (easy  axis):  /3||  =  tt 

1  bias  along  i^  (hard  axis):  =  7r/2 

t  ideal  single-domain  threshold 

c  critical  value  for  dynamic  locking: 

coherent  threshold 


II.  PARAMETERS  OF  A  TYPICAL  PERMALLOY  FILM  (TPF) 


Parameter 

Symbol 

Value 

Gyromagnetic  ratio 

y 

1. 

,87  X 

7 

10  (oe  sec) 

Saturation  magnetization 

M 

o 

7. 

,96  X 

2 

10  gauss 

Exchange  constant 

A 

1( 

3  e 

-1 

rg  cm 

Uniform  uniaxial  anisotropy 

K 

o 

1. 

3  -3 

,6  X  10  erg  cm 

Root-mean-square  local  anisotropy 

5 

4  -3 

X  10  erg  cm 

Mean  crystallite  boundary  spacing 

r 

o 

1, 

.25  X 

10  cm 

Half-thickness 

L 

5 

X  10' 

-  6 

cm 
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APPENDIX  B 

COUPLING  OF  RIPPLE  COMPONENTS 


Here  we  carry  the  iterative  procedure  that  was  used  to  find  the  equilibrium  state  of  the 

magnetization  one  step  further.  We  shall  find  that  the  correction  thus  introduced,  which  arises 

from  coupling  of  the  ripple  components  <p_^  {k  ^  0)  by  random  anisotropy  forces  (but  not  from 

k 

contributions  to  the  torque  nonlinear  in  <p which  are  still  assumed  negligible),  is  small 

k 

provided  fluctuations  of  the  magnetization  from  its  mean  direction  are  small  or,  more  precisely, 
provided 


k 


(B-1) 


k^r^O 


of  the  form  f  f_^  where  f  is  the  Fourier  transform  of  a  random  function  F{r)  satisfying 
k  k*  -k  k 


{III-9a-c).  To  find  this  we  first  obtain  a  generalization  of  the  autocorrelation  function  [see 
(III-IO)] 


k 


{B-2) 


Following  the  same  argument  used  to  derive  c^(r),  we  find 


k,  0 


=  c  (r)  6 
o 


k,  0 


(B-3) 


and  using  the  inverse  Fourier  transform 


k 


{B-4) 


we  have  the  result 


It  k' 


k  k;o 


(B-5) 
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We  now  turn  to  (III-23b)  for  the  ripple  components  <p_^,  which  we  write  in  the  form 

k 


(n+1)  k  1 

■  2U_ 

k  k 


k-k'  k 

k'^0 


(n) 


(k  ^  0) 


(K-6) 


where 


U  =  U  =  K  cos  Zip  +  M  H  cos  {B  ~  (p  )  +  Ak^  +  ZttM  ^  x  sin^(^  —  (p  ) 
k  °  ^  o  ^o  o  ^o 


K  A(^  +  gj 

^  k 


(B-7) 


and  is  the  n^^  approximation  to  <p_^.  Starting  with 


„<0)=0 


{B-8) 


we  rederive  (III-35): 


(1) 


2U 


(B-9) 


Then  the  next  iteration^  gives 
J2)  _  1  / 


2U 


p-.-  Z 


q _ > 

k-k'  k' 


k^ 


k'^0 


(B-10) 


with  mean-square  amplitude 


|^(2)|2^  |^(1)|2^^ 

k  k  4U' 


z 


k -k'  k'  -k 


U 


+  c.c. 


k  *“^k'^0 


q-H.  -.q^  -p^p  - 
k-k'  k"-k  k'  -k' 


k'  k" 


(B-11) 


To  find  the  ensemble  averages  in  (B-11)  we  first  observe,  with  the  help  of  (B-4),  that 


q_^  P_^  P  be  written  as  a  sum  of  terms  each  of  which  has  a  factor 

k-k'  k'  -k 


Q  P  P 
"^1  "^2  "^3 


{B-12) 


t  At  this  level  of  opproximation  there  is  olso  o  correction  to  (p  ,  through  the  term  —2  I  q  (p  of  (lll-23o), 
.  (1)  .  .  (2) 

where  here  is  (p  .  This  leods  to  o  correction  to  (p  which,  however,  may  be  neglected  since  it  is  non- 


lineor  in  (p 


0) 
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where  =  Q(r )  with  r  in  cell  etc.  However,  from  (III-5a,  b)  we  see  that 


P  =  Q  =  QP  =  QP^  =  0  (B-13) 

and  therefore  each  factor  (B-12)  vanishes  (whether  or  not  any  of  the  are  equal)  so  that 


CU  _  P_^  P  _.  =  0 
k-k’  k’  -k 


(B-14) 


Similarly, 


Q  Q  P  P  =  Q  P  6  6 

^3  ^4  ^3' ^4 


(B-15) 


so  that  the  last  ensemble  average  in  (B-11)  may  be  split: 


cu  _  q_.  _p_  p_  =  q^_q_.  _p_p_. 

k -k’  k'’-k  k’  -’k*'  k—k’  k’'~k  k’  — k’ 


I  _  I  ^  I P-.  I  ^  ^ 

k  -k'  k'  k",k' 


(B-16) 


where  we  have  used  (B-5)  for  the  second  step  in  (B-16).  Inserting  (B-14)  and  (B-16)  in  (B-11), 
we  have  the  result 


_  _  leu  _  l"  Ip-.  I" 


(B-17) 


k  k' 

But  now  let  us  assume  (to  be  verified  below)  that  the  main  contribution  to  the  sum  in  (B-17) 
is  from  components  with 

-1 


k’  «  r 


(B-18) 


2  2  2 

Then  we  see  from  (III-14)  that  |  q_^  _^  |  is  independent  of  k',  and  since  P  =  Q  ,  (B-17)  may  be 

k-k' 

written 


-  -  I P-.  I  ^  I  p_  I  ^ 

i„y)|2=  |^H)|2,^^  E  ^ 

k  k  4uE  j.  uj 

k  ^  k’ 


(B-19) 


Finally,  using  (B-9)  and  (B-1),  we  obtain  the  main  result  of  this  appendix 


(B-20) 


where  6  is  the  magnetization 
justified  for  exactly  the  same 


dispersion  with  (p_^  =  .  The  assumption  (B-18)  is  then 

k  ^2 

reasons  that  the  term  (r^k)  was  neglected  in  (III-60). 
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We  conclude  that 


<P 


(2)|2  ^ 

k 


(1)|2 


(B-21) 


and  it  seems  unlikely  that  further  iterations  will  change  |  (p^  |  by  very  much.  Thus,  unless 

k 

(B-1)  is  violated  (in  which  case  the  entire  theory  breaks  down  since  the  nonlinear  torque  terms 
will  no  longer  be  small),  the  assumption  we  have  made  —  that  ripple  components  are  uncoupled 
in  static  equilibrium  (and  also  dynamically  —  see  p.  43)  —  is  indeed  correct. 
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APPENDIX  C 
ON  THE  WKB  METHOD 


Consider  the  general  linear,  homogeneous,  second-order  differential  equation,  which  may 
be  put  in  the  normal  form 

+  F(x)  (^  =  0  .  (C-1) 

We  first  show  that  a  pair  of  solutions  to  (C-1)  is 

<p{x)  =  h(x)  J^|^[g(x)]  (C-2) 


provided  a  generating  function  g(x)  can  be  found  for  the  potential  F(x)  such  that 

c’/  >  /  i%2  3  ,g",2  ,  1  g'"  ,  ,1  2,  ,g\2 

F(x)  =  (g')  -  ^  (^)  +  -  ^  +  (-  -  p  )  (^) 

/  .v2  r" 

=  (g')  - 


where 


h  is  given  by 


l-2r^l/2 


(1  V 


h(x)  =  ■ 

The  proof  requires  only  the  Bessel  function  recursion  relations 

—  J  (g)  =  J  Ag)  +  J  ^Ag) 


2  4-  J  (g)  =  J  <(g)  -  J  ^-(g)  ; 

differentiating  (C-2)  twice  with  the  help  of  (C-5a,b),  we  find 

2n 


h"  -h(g')^  +  [h'g'  + 1  hg"  - (g')^i 


SO  that  (C-1)  is  satisfied  if 

2n 


h”  -h(g')‘ 


'  g  ' 


-Fh 


h'g'  +  I  hg"  -  (g')^  =  0 


(C-3a) 

(C-3b) 


(C-3c) 


(C-4) 


(C-5a) 


(C-5b) 


(C-6) 


(C-7a) 


(C-7b) 


Equation  (C-7b)  may  be  written 


(C-8) 


which  we  integrate,  obtaining  (C-4)  (where  a  multiplicative  constant,  clearly  of  no  interest, 
has  been  set  equal  to  unity).  Finally,  eliminating  h  from  (C-4)  and  (C-7a)  we  arrive  at  (C-3a), 
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Equation  (C-3a)  has  the  significance  that  every  g(x)  generates  an  exactly  soluble  potential  F(x), 


which  depends  on  r.  If  no  g  can  be  found  for 

a  g  such  that  F  approximates  F  (and  has  the 

a  given  F,  it  may  be  possible  with  some  v  to  find 
same  zeros  and  poles)  in  some  range  of  x.  Then 

^(x)  =  (|j  J±^,[g(x)] 

(C-9) 

will  be  an  approximate  solution  of  (C-1)  in  iR. 
and  their  potentials: 

We 

give  a  few  examples  of  generating  functions 

(a)  g{x)  =  Cx^ 

(k  4  0) 

F(x)'  =  (Ctix^"^)^ 

for 

1 

■  2p 

(b)  g(x)  =  C 

F{x)  =  KL^{C^  -  v^) 

(c)  g{x)  =  C  Inx 

2  1-2 

F(x)  =  (C'"  +  i) 

for 

1 

"  =  2 

=  (|-  -  (x  In  x)"^ 

for 

We  also  note  the  asymptotic  behavior  of  the  solution  (C-2),  which  follows  from  the  asymptotic 
expansion  of  the  Bessel  functions 


g  =  I  g I  (C-lOa) 

(c-iob) 

we  find 


<^(x)  -*  [g  “  -TT  <  0  <  TT  (C-lla) 

exp  [7ri(±i^  +  ^)]  cos  [g  +  -|  +  -I)  ]  0  <  0  <  Ztt  (C-llb) 

exp  [27ri(±i^  -i-  y)]  cos  [g  —  y  +  y)]  tt  <  0  <  3;:  ,  (C-1  Ic) 

V  TT  g  ^  L  L 

etc. 

This  approach  to  the  solution  of  (C-1)  is  especially  appropriate  if  the  potential  is  of  the 
WKB  form 

F(x)  =  [Af(x)]^  (C-12a) 

where 


A  »  1  (C-12b) 

and  f(x)  and  its  derivatives  are  0(1)  in  some  range  of  x  away  from  "free"  points  f(x)  =  0. 

In  particular,  if  f(x)  is  slowly  varying  near  a  free  point  or  has  other  unusual  behavior  there. 
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i.e.,  an  inflection  point,  the  standard  methods  may  fail,  and  the  procedure  we  now  give  may  be 
useful.  For  x  in  5R  we  assume  there  exists  a  solution  of  the  form 

oo 


<p{x)  =  exp 


(C-13) 


—  1  —  2 

where  x^  is  a  constant.  Letting  p(x)  =  Pq(x)  +  \  p^(x)  +  .  .  . ,  and  dropping  terms  of  0(X  ),  we 

substitute  (C-13)  and  (C-12a)  in  (C-1)  to  obtain  the  usual  WKB  asymptotic  solutions 


(p{x)- 


X  in  Vf 


exp 


\  f(^)  di 

X 


(C-14) 


Comparison  of  (C-14)  and  (C-1  la)  shows  that  the  asymptotic  value  of  the  generating  function  for 
a  WKB  potential  (C-1 2a)  is 


g(x) - -  X  P  f(|)  d^  =  g(x) 

X  in  IR  '^x 

oo  o 

which  is  easily  verified  by  substitution  in  (C-3a),  giving 
F(x)  =  [Xf(x)]^  [1  +  0(X'^)] 

Now  inserting  (C-15)  into  (C-2)  and  (C-4),  we  find  that 

<;p(x)— —  J^_^(i)=5(x) 

X  in 

oo 


(C-15) 


(C-16) 


(C-17) 


is  an  asymptotic  solution  of  (C-1)  and  (C-1 2a)  for  any  u. 

The  solution  (C-14)  breaks  down,  of  course,  near  a  free  point  f  =  0.  But  if  a  can  be 

~  2 
found  such  that  g(x)  generates  a  close  approximation  to  the  potential  (Xf(x)]  in  a  continuous 

range  of  x  which  includes  both  a  free  point  and  then  (;p(x)  will  be  a  good  approximation  to 

(p{x)  in  that  range.  For  convenience  we  locate  the  free  point  in  whose  neighborhood  we  seek  a 

solution  at  X  =  0,  and  expand  f(x)  about  this  point: 


f(x)  =  x^  f 

j=l 

where  fx  >  — 1.  If  f(x)  is  analytic  at  the  origin,  then 


fx  =  0 


f.  = 

J 


j'  dxJ 


x=0 


(C-18) 


(C-19a) 

(C-I9b) 


and  if  there  is  a  small  residual  potential  at  x  =  0  (a  ’’nearly  free”  point),  a  j  =  0  term  may  be 
included  in  (C-18)  instead  of  the  j  =  1  term.  It  is  usually  advantageous  to  have  the  lower  limit 
x^  in  (C-15)  coincide  with  the  free  point,  and  with  x^  =  0  the  expansion  of  (C-15)  about  x  =  0  is, 
from  (C-18), 
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g(x)  =  Ax' 


fi+ 1 


Z 

J=1 


f. 


hi  +  j  +  1 


(C-20) 


Now  suppose  that  for  a  continuous  range  of  x,  which  includes  both  the  neighborhood  of  the 
free  point  and  part  of  g(x)  may  be  approximated  by  the  j  =  m  term  of  (C-ZO) 


Af 


g(x) 


p  +  m  +  1 


p+m+ 1 


Then  g(x)  generates  a  close  approximation  to  (Af(x)]^  in  !R^  providing 


(C-21) 


1 _ 

2(p  +  m  +  1) 


(C-22) 


as  may  be  seen  from  example  (a),  p.78  .  The  requirement  that  iR^  and  iR^  overlap  is  satisfied 

if  there  is  some  x  in  iR  for  which  condition  (C-lOb)  holds; 
m  m 


2Xl^  |f  I  |x  I  »  y  It  - 

'  m '  '  m '  2  '  4  ' 


(C-23) 


It  may  be  that  no  suitable  JR^  exists;  for  example,  two  terms  of  (C-20)  might  be  required  to 
connect  the  free  point  with  Nevertheless,  it  may  be  possible,  using  (C-3a),  to  find  a  g(x) 

which  approaches  g(x)  for  x  in  IR^  and  generates  a  good  approximation  to  [Af(x)]^  around  the 
free  point. 

To  summarize,  </?(x)  as  given  by  (C-17)  is  an  approximate  solution  to  (C-1)  and  (C-12a)  in 

iR^  and  !R^,  with  g(x)  given  by  (C-15)  with  x^  =  0  and  by  (C-22).  The  asymptotic  formulas 

(C-11)  may  then  be  used  to  establish  the  connection  between  the  coefficients  of  the  solutions 

(C-14)  in  !R^  on  either  side  of  the  free  point.  The  approximate  potential  generated  by  g(x)  differs 

from  (Af(x)]^  by  an  amount  given  by  the  second  term  of  (C-3b),  and  this  may  be  used  to  improve 

57 

the  approximation  if  necessary. 
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APPENDIX  D 

THE  CRITICAL  FIELD 


In  Sec.  VI  we  obtained  a  critical  field,  H  =  H,.h  ,  which  we  identified  as  the  theoretical 

pc  K  pc 

pulse  threshold  for  coherent,  high-speed  switching  Equations  (VI-31)  and  (VI-29b,  c),  which 

determine  h  ,  cannot  in  general  be  solved  in  closed  form.  In  this  appendix  we  use  graphical 
pc 

and  numerical  methods  to  find  h^^  for  the  two  field  configurations  of  greatest  interest. 


1.  HARD-AXIS  BIAS  FIELD  -  EASY -AXIS  PULSE  FIELD 


In  this,  the  most  frequently  used  arrangement, 


~  arcsin  hj^ 

and  therefore  from  (VI-29b,  c)  and  (VI-23-1) 

X  =  h||  -  f  R 

||c  X 

Y  =  h,  -  f  R 

i  y 


where  the  reaction  torque  coefficient  R  is  given  by  (VI-22c)  and 
,  ,1/3,  4/3  -5/6 

X  ||c  i  ^1 

f  =  h  I  /  h  /  g  ' 
y  ||c  1 

gl  -*^1  =  coscp^iO) 


A  necessary  condition  for  Eq.  (VI-31)  to  determine  hi 


is  (p  >  <p  (0),  or 
^oc  ^o 


Y  >  h 


3 

i 


(D-la) 

(D-lb) 

(D-lc) 

(D-2a) 

(D-2b) 

(D-3a) 

(D-3b) 

(D-3c) 


which  leads  to 


R  <  h 


2/3 

i 


gV2 


where  we  have  used  the  relation 


f 

X 


=  tan  (D  (0) 
o 


(D-4) 


(D-5) 


We  shall  assume  that  (D-4)  is  also  sufficient  for  (VI-31)  to  be  valid,  which  is  equivalent  to 
assuming  that  «  1.  Equation  (VI-31)  is  most  conveniently  solved  for  R  (given  and  hj^) 
by  the  graphical  method  sketched  in  Fig.  D-1,  where  we  have  made  use  of  (D-5).  In  Fig.  D-2 
we  plot  vs  hj^  with  R  as  a  parameter. 
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Fig.  D-1 .  Graphical  canstructian  used  to  find 
reaction  torque  coefficient  R  as  a  function  af 
hord-oxis  bias  field  and  critical  easy-axis 


switching  field  h 


Ih 


Fig.  D-2. 
theory  and 


High-speed  switching  threshold  curves  with  reaction  torque  coefficient  R  as  a  parameter: 

experiment  (h  ^  =  h  , ,  with  hard-axis  bias  field  h , ). 

p3  1 1 c  1 
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For 


P  .  .  2/3  3/2 

h||^  is  given  in  first  approximation  by  (VI-33),  which  becomes 


R 


3/2 


c  .  5/4 


(D-6) 


(D-7) 


However,  this  expression  may  be  seriously  in  error  since  it  neglects  the  growth  time  of  the 
reaction  torque;  we  show  it  as  a  dashed  line  in  Fig.  D-2. 

II.  ZERO  BIAS  FIELD 

In  this  arrangement  (which  is  more  frequently  used  to  study  quasistatic  than  pulse  switching) 

(D-8a) 
(D-8b) 


-  TT  -/s  {0  <  (3^ 


h  -*  h 
P 

<^^(0)  =  0 


h,  =  0 
b 


and  therefore 


X  =  h  cosiS 
c  ^ 


Y  =  sin/?  —  (h^  sin/?)  '  R 


The  condition  for  (VI-31)  to  determine  h  is  now 

c 


R  <  (tanjS) 


2/3 


(D-9a) 

(D-9b) 

(D-10) 


We  solve  (VI-31)  for  h^  cos /3  as  a  function  of  h^  sin^S  and  R,  and  the  result  (with  R  a  parameter) 
is  plotted  in  Fig.  D-3. 

For 


R  >  (tan 3’)^'^^ 

(VI-33)  gives  the  first  approximation 
,3/2 


h  = 
c 


R" 


(D-11) 


(D-12) 


sin)3 

which  is  shown  as  a  dashed  line  in  Fig.  D-3. 
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Fig.  D-3.  High-speed  switching  threshold  curves  with  reoction 
torque  coefficient  R  os  porometer  (hp3  =  h^  at  an  ongle  P  to 
eosy  direction  opposite  initial  magnetizotion). 
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